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CHAPTER  I 


1.0  INTRODUCTION 

To  demonstrate  Maxwell's  theory,  that  light  waves  and  radio  waves  will  behave  is 
a  similar  fashion,  Heinrich  Hertz,  in  1886,  experimentally  reflected  radio  waves  (66cm)  off 
of  a  metallic  surface.  By  demonstrating  that  light  waves  and  radio  waves  are  both  part  of 
the  electromagnetic  spectrum  and  can  be  reflected  off  metallic  or  dielectric  bodies,  Hertz 
set  the  stage  for  radar  technology[ll.  The  first  RADAR  (RAdio  detection  And  Ranging) 
systems  were  used  to  warn  of  approaching  hostile  aircraft  in  the  1930’s  but  with  the 
invention  of  the  laser  in  1960,  distance  and  velocity  could  be  measured  more  accurately 
because  of  the  shorter  wavelengths  that  light  provides  over  those  provided  with 
conventional  radar. 

The  first  laser  radar  systems  were  know  as  LIDAR  (Light  detection  And  Ranging) 
systems.  Although  it  was  recognized  that  the  Light  was  laser  light,  work  has  been  done 
with  Xenon  and  other  flash  lamps  so  the  name  has  been  changed  to  LADAR  (LAser 
detection  And  Ranging)  [2]  when  referring  specifically  to  laser  light.  Laser  range  finders 
were  first  used  in  many  aspects  of  the  Vietnam  war;  ground  troops  used  hand  held  lidars, 
lidars  were  mounted  on  tanks  to  aid  in  measuring  distances,  and  "smart"  bombs  used  lasers 
to  track  targets.  Other  applications  for  laser  radars  include  remote  target  tracking,  aircraft 
altimeters,  and  atmospheric  studies  including  wind  profiling  and  pollution  monitoring  [3]. 
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Ladar  wavelengths  commonly  used  today  are  10.6|im  (CO2),  1.06|im  (Nd:YAG), 
and,  recently,  2.09|i,m  (Tm,  Ho:YAG),  where  the  information  in  parentheses  indicates  the 
appropriate  gain  medium  used  for  the  corresponding  wavelengths.  In  196X,  Raytheon 
demonstrated  the  first  coherent  detection  ladar  system  [4]  using  a  CO2  laser.  Coherent 
detection  differs  from  direct  detection  in  that  the  return  iiequency  shifted  signal  is  mixed 
with  a  reference  signal  called  a  local  oscillator.  This  produces  a  beat  signal  that  is  easier 
to  detect  than  just  the  return  from  a  target.  CO2  ladars  are  used  in  the  coherent  mode  of 
operation  because  of  the  background  noise  produced  by  the  environment  at  10.6  pm.  To 
partially  reduce  this  background  noise,  the  optics  used  in  a  CO2  .sy.stem,  along  with  the 
detectors,  have  to  be  cooled  with  liquid  nitrogen.  This,  along  with  the  heavy  rf  power 
supplies,  make  CO2  systems  heavy  and  expensive  to  operate  [5]. 

Nd:YAG  ladar  systems  are  being  used  because  they  are  solid  state  sy.stems.  They 
do  not,  however,  have  the  efficiencies  or  the  output  energies  of  the  10.6  |im  CO2  systems, 
though,  as  they  are  solid  state  systems  they  are  light-weight  and  less  expensive  to  produce. 
There  are  also  efficient  detectors  at  this  wavelength  that  do  not  require  cooling  with  liquid 
nitrogen  and  the  laser  source  is  continuously  tunable  over  a  few  nanometer  region  [5]. 
CO2  systems  are  also  tunable  but  only  to  other  discrete  molecular  lines.  All  solid  state 
1.06pm  ladar  systems  were  direct  detection  systems  until  1986  when  Stanford  University 
.successfully  demonstrated  coherent  detection  of  signals  from  clouds  and  atmo.spheric 
aerosols  particles  [6],  A  more  efficient  second  generation  Nd:YAG  ladar  system  was  built 
in  1988-89  by  Coherent  Technologies,  Inc.  (CTI)  in  Boulder,  Colorado.  The  experimental 
design  is  such  that  the  system  is  more  operational  and  the  data  analysis  is  more  real  time 
than  the  system  built  by  Stanford  University.  CTI  has  used  their  system  to  do  atmospheric 
absorption  studies  and  wind  profiling.  At  this  ti.me  it  is  the  only  known,  mobile  coherent 
1.06  pm  Nd:YAG  ladar  system  [7]. 

In  1987  the  U.S.  Air  Force  approached  CTI  about  designing  an  eye  safe,  coherent, 
solid  state  ladar  system.  The  primary  desire  for  an  eye  safe  system  is  so  that  a  pilot  will 
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not  blind  his  wing  man,  or  other  personnel,  when  he  has  his  ladar  system  engaged.  This  is 
just  an  example  of  a  scenario  in  which  it  would  be  desirable  to  have  an  eye  safe  system. 
The  wavelength  chosen  by  CTI  was  2.09  |im  because  materials  have  been  developed  that 
emit  in  the  2.09  |xm  region.  These  materials  can  be  efficiently  pumped  using  laser  diodes, 
and  there  are  high  quantum  efficiency  detectors  as  well  as  a  high  atmospheric  tran.smission 
in  that  region  [6].  Along  with  CTTs  1.06  |im  system  their  2.09  pm  system  is  u.sed  to  do 
wind  profiling,  range-resolved  wind  velocity  and  aerosol  backscatter  measurements.  Until 
now  CTI  had  the  only  operating  2.09pm  coherent  ladar  system  [6]. 

CTI  then  released  its  research  to  Wright  Research  and  Development  Center  at 
Wright-Patterson  Air  Force  Base,  Ohio  in  1990.  Their  2.09  pm  system  has  been 
duplicated  and  is  being  used  to  compare  coherent  (i.e.  heterodyne)  detection  to  incoherent 
(i.e.  direct)  detection  at  2.09pm.  In  1966,  J.  W.  Goodman  compared  the  detection 
techniques  for  wavelengths  shorter  than  1pm.  He  found  that  when  there  is  a  large  amount 
of  background  noise,  coherent  detection  systems  have  greater  sensitivity.  Coherent 
systems  also  perform  better  than  incoherent  systems  when  high  velocity  resolution  in 
needed.  However  incoherent  systems  perform  with  higher  signal-to-noise  ratios  than 
coherent  systems  when  operated  in  low  noise  environments,  such  as  space  based 
applications,  and  at  small  rates  of  false  alarm  [8].  Robert  J.  Keyes'  analysis  in  1986  found 
that  incoherent  detection  could  perform  with  nearly  the  same  signal-to-noise  ratio  or  even 
higher  signal-to-noise  ratio  than  coherent  detection  for  wavelengths  in  the  near-infrared 
and  visible  regions.  In  the  mid-  to  far-infrared  regions  the  detection  schemes  gave  similar 
results  [9].  A  comparison  of  the  detection  techniques,  however,  has  never  been  done  for  a 
coherent  solid  state  system  and  it  has  never  been  done  specifically  at  2.09pm  which  is  in 
the  mid-infi'ared  region. 

The  comparison  of  the  detection  techniques  will  be  done  using  a  2pm  ladar  system 
that  is  capable  of  both  coherent  and  incoherent  modes  of  operation.  The  comparison  Is 
made  based  on  the  probability  of  detection  for  each  detection  scheme  as  an  analysis  of  this 
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type  has  not  yet  been  performed  for  an  eye  safe  coherent  solid  state  ladar  system.  The 
probability  of  detection  is  a  function  of  the  statistics  of  detected  light  returning  from  a 
distant  target.  These  statistics  are  dependent  on  the  type  of  target,  the  detection  scheme, 
the  electronics  in  the  detection  circuit  and  the  types  of  noise  in  the  system.  An  analysis  of 
these  factors  is  pursued  in  this  text  to  compare  tlie  detection  teclmiques. 


CHAPTER  II 


2.0  Eve  safe  LADAR  Rfeauirements 

For  a  system  to  be  considered  eye  safe,  direct  exposure  to  the  transmitted  laser 
beann  must  not  damage  the  eye  under  normal  conditions.  Different  parts  of  the  eye  are 
sensitive  to  different  wavelengths  of  light  For  example,  the  retina  is  sensitive  to  visible 
light  (400-700nm)  and  ER-A  (700-1400nm)  radiation,  whereas,  the  lens,  aqueous  humor 
and  cornea  absorb  UV  (200-400nm),  IR-B  (1.4-3|i.m)  and  IR-C  (3-1000|i,m)  wavelengths. 
At  2.09|xm,  the  cornea  absorbs  about  75%  of  the  incident  energy,  while  the  remaining 
25%  is  absorbed  by  the  aqueous  humor.  The  primary  mechanisn-is,  then,  by  which  the  eye 
may  be  damaged  by  2.09|im  radiation  are  excess  heat  generation  in  the  mostly  water  based 
aqueous  humor  and,  more  importantly,  the  formation  of  comeal  cataracts  [8].  Limits  must 
therefore  be  set  with  regard  to  exposure  duration  and  intensity  so  as  to  minimize  eye 
damaging  effects. 

The  American  National  Standards  Institute  (ANSI)  has  issued  standards  for 
maximum  permissible  exposure  (MPE),  which  is  defined  as  the  radiant  exposure  which 
individuals  may  receive  without  harmful  biological  effects  [8].  Since  damage  to  the  eye 
depends  on  the  wavelength  and  the  exposure  duration,  the  standards  vary  according  to 
wavelength  and  exposure  time.  For  IR-B&C  wavelengths,  the  standards  are  given  in 
Table  1. 

For  our  system  (see  Fig.  9,  which  will  be  more  fiiliy  described  later)  the  beam 
internal  to  the  system  will  damage  the  eye.  Once  the  system  is  aligned  the  only  beam  that 
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Table  1: 


Exposure  limits  for  direct  exposure  to  a  laser  beam  foi  IR-B&C  wavelengths. 


Wavelength 

Exposure  Time 
(t)  seconds 

Exposure  Limits 

10-9  to  10-7 

10"'  y  : 

/cm* 

10-7  to  10 

MlltEIRuSiESRSIlllll 

10  to  3x!04 

0.1'%.. 

escapes  the  system  is  the  expanded  transmitted  beam.  The  exposure  limit  for  this 
transmitted  pulse  of  length  500  ns  is  [8] 


£i,™«)=o.56.V7x.. 

(0=0.56.^500x10^;^,,  ^2.1) 

=  1.489x10-'/, 

/cm 

=  14.9'^,. 

/cm* 

For  our  system  with  a  typical  transmitted  energy  of  15  mJ,  with  an  expanded  /:  beam 
diameter  of  4  cm  at  the  exit  aperture  of  the  telescope,  the  radiant  exposure  an  individual 
would  receive  by  glancing  into  the  exit  aperture  of  the  telescope  can  be  calculated  as 


(Output  Energy) 
(  Area  of  Beam) 


\5mJ 

4^ncm^ 


=  0.298'^/„,. 


(2.2) 


for  a  single  pulse.  This  is  only  2%  of  the  single  pulse  MPE  of  14.9mJ/cm2  for  this  system, 
according  to  ANSI  standards  [8]. 

In  order  to  assess  the  effects  of  extended  exposure  to  2.09|.im  radiation,  the  la.ser 
hazard  assessment  program  LHAZ,  developed  by  Armstrong  Laboratory  according  to 
ANSI  standard  136.1-1986,  has  been  used  [9].  According  to  this  program,  an  individual 
could  stare  directly  into  the  e.xit  aperture  of  the  telescope  of  our  system  (see  Figure  9)  for 
8.3  hours  and  only  receive  68%  of  the  extended  MPE  of  7 1 .5  J/cm^.  Ba.sed  then  on  .single 
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pulse  and  extended  exposure  limits  our  2.09|im  LADAR  system  is  considered  to  be  eye 
safe  with  respect  to  the  transmitted  laser  radiation. 

By  contrast,  using  the  same  pulse  duration,  beam  diameter  and  pulse  energy,  the 
single  pulse  MPE  is  given  by  ANSI  to  be  5  pJ/cm^  for  1.06  pm  radiation  [9].  The  actual 
radiant  exposure  from  a  single  pulse  (0.298  mJ/cm^)  would  thus  be  enough  to  damage  the 
eye  at  this  wavelength.  For  this  reason,  a  comparable  1.06  pm  LADAR  system  would  not 
be  considered  eye  safe.  However,  the  transmitted  beams  of  CO2  LADARs,  under  the 
same  conditions  as  considered  for  the  2.09  pm  system,  are  eye  safe  since  the  single  pulse 
MPE  is  the  same  as  for  2.09  pm  [9],  as  can  be  seen  in  Table  1.  As  discussed  previously, 
though,  CO2  systems  do  not  have  the  tuneability  of  the  2.09  pm  systems  allowing  for 
tuning  around  atmospheric  absorption  lines  and  unlike  solid  state  systems  they  have  to  be 
cooled  with  liquid  nitrogen. 


CHAPTER  in 


a.O  THEORY 

The  comparison  of  coherent  to  incoherent  detection  will  be  performed  assuming 
both  a  speckle  and  a  glint  target.  The  basis  for  this  comparison  will  be  the  probability  of 
detection.  To  understand  the  probability  of  detection,  the  role  of  the  statistical  nature  of 
the  noise  and  the  return  signal  plus  noise  must  first  be  understood.  The  noise  distribution 
shown  in  Fig.  1  [1]  is  the  probability  density  function  (PDF)  that  represents,  generally,  the 
fluctuation  of  the  noise  current  in  a  radar  system  when  there  is  no  target  present. 


Figure  1:  Probability  density  functions  representing  both  signal  and  noise  distributions, 
with  graphical  representation  of  the  probabilities  of  detection,  false  alarm  and  miss. 
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The  shape  and  the  position  of  this  probability^  distribution  is  dependent  on  the  detection 
technique  used  and  on  the  average  value  of  the  noise  current,  /yvo/jre-  When  a  signal  is 
present,  the  average  value  of  the  signal  current  Isignal  generally  greater  than  the 
average  value  of  the  noise  current,  so  that  the  PDF  representing  the  fluctuation  in  the 
combined  return  signal,  plus  noise  current,  is  centered  about  an  average  value  of  the  return 
signal  plus  noise.  The  second  distribution  shown  in  Fig.  1  represents,  generally,  the  sum 
of  the  signal  and  noise  currents. 

In  Older  to  decide  whether  a  value  measured  by  the  detector  is  from  an  actual 
target  or  whether  it  is  noise,  a  threshold  current  fThresholcl  shown  in  Fig.  1. 

Wlienever  a  current  produced  by  the  detector  is  larger  than  the  threshold  current,  a  target 
is  said  to  have  been  detected.  Whether  or  not  a  real  target  exists  has  yet  to  be  detennined. 
The  probability  that  a  target  has  been  detected  is  called  the  probability  of  detection,  which 
is,  mathematically,  the  area  under  the  PDF  for  the  signal  and  noise  greater  than  the  set 
threshold.  There  is  also  the  probability  that  the  current  produced  by  the  detector  exceeds 
the  threshold  due  to  noise  effects  only.  The  probability  of  this  occurring  is  called  the 
probability  of  false  alarm  and  is,  mathematically,  the  area  under  the  noise  current 
distribution  greater  than  the  set  threshold  current.  And  finally,  tliere  is  a  small,  undesirable 
possibility  that  a  signal  from  a  real  target  can  produce  a  detected  current  less  than  the 
threshold  current,  causing  the  target  not  to  be  detected.  This  is  called  the  probability  of 
miss. 

In  order  to  find  the  probability  of  detection  for  the  2.09  pm  LADAR  system  the 
average  value  of  the  dominating  noise  source(s)  needs  to  be  determined  as  will  be 
discussed  in  the  following  section.  The  PDF  that  the  dominating  noise  source(s)  exhibit 
will  then  be  discussed  for  both  coherent  and  incoherent  detection.  The  probability 
distributions  of  the  noise  will  then  be  used  in  conjunction  with  desired  probabilities  of  false 
alarm  to  find  the  needed  threshold  for  that  false  alarm  rate.  The  PDF  representing  the 
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fluctuations  in  the  detector  output  caused  by  the  type  of  target  and  the  detection  process 
will  then  be  identified  for  a  glint  and  a  speckle  target  for  both  coherent  and  incoherent 
detection.  Once  these  four  distributions  are  identified,  the  threshold  values  found  when 
calculating  the  probabilities  of  false  alarm  will  be  used  with  the  PDFs  representing  the 
fluctuation  in  the  detector  output  current  to  calculate  the  corresponding  probabilities  of 
detection. 


3.1  Noise  Sources 

When  trying  to  detect  a  signal  there  is  always  noise  present.  Electrically,  noise  is 
expressed  as  a  mean  squared  current  or  voltage  fluctuation  around  a  DC  value,  which  Is 
also  called  the  variance.  This  noise  may  be  generated  by  the  randomness  associated  with 
the  detection  of  a  signal  (shot  noise),  the  leakage  current  from  the  detector  (dark  current 
noise),  or  the  current  flowing  through  a  resistor  in  the  detection  circuit  (thermal  noise  or 
Johnson  noise).  Other  sources  of  noise  are  stray  light  striking  the  detector  that  is  not  from 
the  intended  source  (background  noise),  and  noise  in  the  post-detection  circuit  caused  by 
an  amplifier  (amplifier  noise).  The  following  is  a  discussion  of  these  noise  sources. 

3.1.1  Shot  Noise 

Shot  noise  is  the  fluctuation  in  the  current  due  to  the  discrete  nature  in  which 
charge  carriers  are  produced.  This  fluctuation  can  be  seen  if  a  DC  current  is  looked  at  on 
a  short  time  scale.  This  fluctuation  has  a  Poisson  probability  density  function  where  the 
variance  is  equal  to  the  mean  [10].  The  mean  squared  current  fluctuation  seen  at  the 
output  of  an  electronic  filter  coupled  with  the  detector  is  (see  Appendix  A)  [10] 


(4>  =  2?/B, 


(3.1.1) 
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where  q  is  the  charge  on  an  electron,  /  is  the  average  ciurent  and  B  is  the  electrical 
bandwidth  of  the  detection  circuit  determined  by  the  electrical  filter, 

3.1.2  Dark  Current  Noise 

The  dark  current  of  a  detector  is  the  leakage  current  produced  when  there  is  no 
energy  incident  of  the  surface  of  the  detector.  TTtis  leakage  current  is  always  there  and  the 
value  of  the  leakage  current  differs  from  detector  to  detector  even  if  the  detectors  were 
produced  in  the  same  batch.  The  equation  for  the  mean  square  value  of  the  dark  current 
noise  is  the  same  as  the  equation  for  shot  noise  [11]  except  is  used  to  represent  the 
average  dark  current ,  that  is 

{ll,)  =  2ql^,B  (3.1.2) 


3.1.3  Johnson  Noise 

Johnson  noise  or  thermal  noise  is  fluctuation  caused  by  the  thermal  motion  of 
charge  carriers  in  a  dissipative  element  An  example  of  a  dissipative  element  that  produces 
thermal  noise  is  the  resistive  load  used  to  measure  the  signal.  The  equation  used  to 
calculate  the  mean  squared  Johnson  noise  is  [1 1] 


AkTB 

^Load 


(3.1.3) 


where  k  is  Boltzman's  constant,  T  is  the  temperature  of  the  element  in  Kelvin,  and  R^oad 
is  the  resistance  of  the  load  as  seen  by  the  detector. 
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3.1.4  Background  Noise 

Another  source  of  noise  in  a  ladar  system  is  background  noise.  Background  noise 
is  the  shot  noise  produced  by  the  detection  of  radiation  that  has  reflected  off  or  comes 
from  the  earth,  clouds,  the  atmosphere  and  the  sun.  The  shot  noise  current  produced  by 
the  background  is  [4] 

{il)  =  2qBP,,R^„.  (3.1.4) 

where  is  the  responsivity  of  the  detector  and  is  the  power  incident  on  the 
detector  produced  by  solar  backscatter.  The  power  from  solar  backscatter,  P^g,  has  been 
shown  to  be  [4] 

^sB  ~  ^R  ’  (3.1.5) 

where  ki  'is  the  fraction  of  the  solar  radiation  that  penetrates  the  Earth's  atmosphere, 
is  the  solar  irradiance,  Qr  is  the  solid  angle  over  which  energy-radiates  from  the  radiating 
body,  p  is  the  target  reflectivity,  11575  is  tlie  optical  efficiency  of  the  system,  is  the  area 
of  the  receiver,  and  AX  is  the  wavelength  band  of  the  optical  bandpass  filter  to  be  placed 
directly  in  fiont  of  the  detector,  which  will  be  centered  around  the  wavelength  of  interest. 
Such  a  filter  will  eliminate  wavelengths  other  than  the  laser  wavelength  of  interest,  thus 
decreasing  the  background  noise. 

3.1.5  Amplifier  Noise 

Amplifier  noise  is  the  noise  added  to  the  signal  through  the  process  of 
ampJification.  The  noise  added  to  a  system  by  an  amplifier  is  either  defined  as  a  noise 
equivalent  temperature  or  a  noise  power  spectral  density.  When  the  noise  is  specified  as  a 
equivalent  noise  temperature  the  noise  current  added  by  the  amplifier  can  be  calculated  by 
using  Eq.  (3.1.3)  which  is  the  equation  for  Johnson  noise.  The  equivalent  noise 
temperature  is  used  as  the  temperature  instead  of  room  temperature  which  is  used  in  most 
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cases.  When  the  noise  power  spectral  density  is  given  it  is  usually  given  in  units  of 
To  find  the  noise  current  added  by  the  amplifier  the  noise  spectral  density  is  s(quured, 
multiplied  by  th?  bandwidth  and  then  divided  by  the  square  of  the  input  impedance.  These 
calculations  are  shown  in  Chapter  5  for  specific  specified  values  given  by  the  manufacture. 


3.2  Probability  Distributions  of  Noise 

The  fluctuations  of  current  in  the  detection  circuit  are  caused  by  the  detection 
process  and  other  noise  sources  as  described  earlier.  These  fluctuations,  being  random, 
can  be  expressed  using  probability  density  functions  (PDF's).  Both  coherent  and 
incoherent  detection  techniques  have  PDFs  representing  the  probability  distribution  of  the 
primary  noise  source.  Using  these  distributions,  equations  for  the  probability  of  false 
alarm  can  be  found. 


3.2.1  Probability  Distribution  of  Noise  for  Incoherent  Detection 

For  incoherent  detection  the  dominating  noise  source,  as  will  be  shown  later  in 
section  5. 1 .5.2,  is  return  signal  shot  noise.  Tlie  shot  noise  current  fluctuation  around  a 
DC  average  current  is  a  Poisson  PDF  but  at  these  high  event  densities  it  is  approximated 
as  a  Gaussian  PDF  with  a  zero  mean  expressed  as  [1,12] 


Pin{i)  = 


(3.2.1) 


where  i  is  the  instantaneous  detector  output  current  and  {if„cokSN)  is  the  mean  squared 
dominating  noise  current,  which  for  this  case  is  signal  shot  noise. 
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3.2.2  Probability  Distribution  of  Noise  for  Coherent  Detection 

For  the  heterodyne  or  coherent  detection  case,  the  desired  doininating  noise  is 
local  oscillator  (LO)  shot  noise.  Shot  noise  can  be  represented  by  a  Gaussian  PDF  before 
it  is  peak  envelope  detected,  which  is  the  common  method  of  analysis  used  when  doing 
coherent  detection  [12].  Once  peak  envelope  detected.  Rice  [14]  has  shown  that  the  noise 
current  envelope  follows  a  Rayleigh  PDF  Pco(i)  given  by  [12] 


(3.2.2) 


where  the  dominating  mean  squared  noise  current  for  coherent  detection  is  LO  shot  noise, 
(ir„h.sN ) » ^  discussed  later. 


3.3  PDF’s  Representing  the  Combination  of  Signal  and  Noise 

Not  unlike  the  noise  in  the  detection  process,  the  return  energy  also  fluctuates  and 
the  combination  of  that  fluctuation  and  the  fluctuation  of  the  noise  in  the  system  can  be 
represented  by  appropriate  probability  density  functions  (PDF's).  These  distributions  will 
be  different  depending  on  the  detection  scheme  used,  the  type  of  target,  and  the  amount  of 
atmospheric  turbulence.  The  targets  of  interest  are  glint  and  speckle  targets. 

This  section  includes  a  description  of  the  different  PDFs  for  each  of  the  detection 
schemes  for  each  type  of  target.  Once  the  PDFs  for  each  case  are  established,  the 
probability  of  detection  can  be  calculated  and  these  detection  techniques  can  be  compared. 
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3.3. 1  The  PDF 's  for  Incoherent  Detection 

For  incoherent  (or  energy)  detection  there  are  different  PDFs  for  the  sum  of  the 
noise  and  the  signal  currents,  for  both  glint  and  speckle  targets.  A  discussion  of  the  PDFs 
for  both  targets  follows. 


3.3.1  A  Glint  Target  returns,  or  specular  reflections,  have  return  energies  that  can  be 
calculated  because  the  characteristics  of  the  particular  target  are  known.  The  process  of 
detecting  the  return  energy  cau.ses  fluctuations  in  the  output  current.  This  fluctuation  can 
be  de.scribed  by  a  Poisson  probability  distribution  given  as  [12] 


F(^)  = 


(3.3.1) 


where  k  is  the  number  of  photoelectrons  emitted  by  the  detector  and  <k>  is  the  average 
number  of  photoelectrons  emitted  by  the  detector.  An  example  of  a  Poisson  distribution  is 
shown  in  Fig.  2  for  only  a  few  incident  photons.  In  our  situation  however,  there  is  a  much 
larger  number  of  photons  striking  the  detector  producing  a  larger  current.  The  central 
limit  theorem  states  that  as  the  number  of  statistically  independent  events  occur  without 
limit,  the  Poisson  probability  distribution  will  tend  toward  a  Gaussian  distribution 
[2,12,15]  given  generally  as 


p{k)  = 


1 

■  ik-mf] 

V™™  2 

ilncr 

(3.3.2) 


which  has  a  mean,  m,  equal  to  the  variance,  a^,  (see  Appendix  A) 

a^=m.  (3.3.3) 


This  can  also  be  shown  graphically  as  in  Fig  3. 


Probability  Density  ProbabHity  Density 
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Since  the  target  of  interest  for  this  case  is  a  glint  target,  the  return  will  be  very 
large  and  the  Poisson  distribution  can  be  approximated  by  a  Gaussian  distribution.  The 


Figure  2:  Poisson  distributions  with  different  mean  values. 


Figure  3:  Comparison  between  Poisson  and  a  Gaussian  distribution. 
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dominant  noise  for  the  incoherent  detection  case  (signal  shot  noise)  also  has  a  Gaussian 
distribution  with  zero  mean  as  stated  earlier.  So  the  combined  distribution  of  the  two 
independent  Gaussian  variables  is  a  Gaussian  distribution  where  the  mean  is  equal  to  the 
sum  of  the  signal  and  noise  current  means  and  the  variance  is  equal  to  the  sum  of  the 
independent  signal  and  noise  current  variances  [16].  The  PDF  representing  the  fluctuation 
in  the  current  i  for  incoherent  detection  with  a  glint  target  Pkj  is  [12] 


Pig(0  = 


1 

■^^‘^{{hncoh.SN  )  )) 

^{{hncoKSN  )  4 

(3.3.4) 


where  is  the  average  current  produced  by  the  return  energy  from  a  glint  target,  q  is 
the  charge  on  an  electron,  B  is  the  bandwidth  of  the  detection  electronics,  and 

is  tlie  mean  squared  signal  current,  (Further  explanation  is  found  in  Appendix  A.) 


3.3  J. 2  Speckle  or  Diffuse  Targets  have  return  energies  that  are  random  due  to  the 
surface  irregularities  in  the  targets.  The  detection  statistics  for  this  case  are  conditionally 
Poisson  conditioned  on  knowing  the  amount  of  return  energy.  Since  the  return  energy  is 
random,  the  unconditional  PDF  that  represents  k  signal  photoelectrons  emitted  is  a 
negative  binomial  distribution  given  by  Goodman  as  [15] 

«• 

=  (3.3.5) 

0 

where  k,  the  number  of  photoelectrons  can  be  converted  to  current  by  multiplying  by 
q/At  where  q  is  the  charge  on  an  electron  and  At  is  the  transmitted  pulse  length  which  will 
be  discussed  later  in  Section  4,1.2,  p(W)  is  the  probability  density  function  of  the  return 
energy  (MO  incident  on  the  detector  during  a  pulse  and  is  given  by  the  Gamma  density  as. 
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a*^tV“~'exp(-alV) 

r(M) 


;1V^0 
;  Otherwise 


(3.3.6) 


where  M  is  the  number  of  spatial  correlation  cells  received  by  the  detector,  a  - 
{w)  is  the  average  return  energy  and  F  is  the  gamma  function.  M  can  be  thought  of  as  a 

measure  of  the  spatial  granularity  of  the  target  as  seen  by  the  receiver.  Furthermore, 
is  conditional  PDF  of  the  total  (signal  and  noise)  photoelectron  count 

emitted  by  the  detector.  The  conditional  PDF  for  the  signal  photoelectron  count  alone 
(Pj(A:|1V))  is  a  Poisson  probability  density  function  given  as 

^s(^F)  = - — - (3.3.7) 


where  T]  is  the  quantum  efficiency  of  the  detector,  h  is  Plank's  constant,  v  is  the  optical 
frequency  and  TiW/Ziv  is  the  mean  number  of  detected  photoelectrons.  Fortunately,  in  the 
presence  of  a  large  photoelectron  count  rate,  a  Poisson  distribution  can  be  approximated 
as  a  Gaussian  distribution.  Thus  the  probability  that  k  photoelectrons  are  emitted  by  the 
detector  is  approximately 


p,(m) = 


1 


/ 

exp 


IbW  J 


(3.3.8) 


where  b  =  T\/hv,  which  is  just  a  constant  chosen  to  simplify  the  equation. 

The  dominant  noise  for  incoherent  detection,  which  is  signal  shot  noise,  also  has  a 
Gaussian  PDF.  The  sum  of  two  Gaussian  random  variables  (signal  and  noise)  has  a 
Gaussian  PDF  where  the  variance  is  equal  to  the  sum  of  the  variances  of  the  random 
variables.  So  the  PDF  representing  the  number  of  electrons  in  the  detector  circuit  is  a 
Gaussian  PDF  of  the  form 


S+N 


(k-{{n)+bW)) 
2{{n)  +  bW) 
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(3.3.9) 


where  (m)  is  the  mean  number  (and  the  variance)  of  noise  photoelectrons.  Now  using 
Eqs.  (3.3.5),  (3.3.6),  and  (3.3.9)  the  distribution  of  the  return  signal  and  noise  is  given  by 


p,E(i‘)=lPs,Ab\w)p(w)dw 


1 

■  {k-{{n)+bW))'‘' 

exp(-aW0 

^2n{{n)  +  bW) 

2{(n)  +  bW) 

r(M) 

dW. 


(3.3.10) 


An  example  of  this  distribution  is  plotted  in  Fig.  5  using  1  as  the  average  noise  ({«)),  2  as 
the  mean  return  signal  ({IV)),  and  4  for  the  number  of  speckle  lobes  (M). 


Figure  4:  This  is  an  example  of  the  distribution  of  the  return  energy  and  the  noise  for 
incoherent  detection  with  a  speckle  target. 


20 


3.3.2  Probability  Distribution  for  Coherent  Detection 

For  coherent  (or  heterodyne)  detection  the  same  probability  density  function  can 
be  used  to  represent  the  signal  current  due  to  either  a  glint  or  a  speckle  target,  or  even  a 
target  containing  both  glint  and  speckle  components.  This  density  is  called  the  Rician 
PDF  and  has  the  form  [2] 

p(i)  =  ^exp 

where  /  represent  the  total  instantaneous  peak  envelope  detected  signal,  is  the  mean 

squared  strength  of  the  speckle  plus  shot  noise  component  expressed  as 
-  {{^LkSN  )  (hiffus,))  *  {^DiffLU  )  represents  the  mean  squared  signal  current  from  a 

diffuse  target,  (4>a,5/v)  is  the  mean  squared  noise  current  (shot  noise)  for  coherent 

detection,  is  the  zeroth  order  modified  Bessel  function  of  the  first  kind,  and  is  the 

portion  of  the  total  signal  arising  from  a  glint  target  component.  An  example  of  this 
combined  distribution  is  shown  in  Fig.  5  where  set  to  5  and  {iliffust)  set 


i^+i^ 

*  ^  '•amt 


(3.3.11) 


2Diffuse  =  Glint 


Figure  5:  A  Rician  distribution  with  different  components  of  diffuse  and  glint  targets. 
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to  10,  while  was  varied  to  give  different  glint  to  noise  ratios.  The  discussion  of  the 
probability  distributions  for  each  target  individually  follows. 


3.32.1  A  SpecklelDiffuse  Target  is  by  definition  optically  rough  and  scatters  incident 
light  randomly.  The  light  scattered  from  each  individual  scattering  centers  interfere  with 
each  other  to  produce  a  speckle  pattern,  which  when  viewed,  resembles  random  light  and 
dark  patches.  The  randomness  in  the  speckle  pattern  and  the  randomness  associated  with 
the  random  phase  of  these  speckles  produces  a  random  fluctuation  in  the  current  produced 
by  the  detector  which  can  be  represented  by  a  Gaussian  distribution  [2]  prior  to  envelope 
detection.  Since  both  the  signal  and  the  noise  currents  are  represented  by  Gaussian  PDF's, 
their  combination  can  be  represented  by  a  Gaussian  PDF  where  the  mean  squared  value  is 
equal  to  the  sum  of  the  mean  squared  values  of  the  signal  and  noise  currents.  With 
coherent  detection,  the  detected  Gaussian  distributed  signal  will  be  passed  through  a 
bandpass  filter  centered  around  the  intermediate  frequency.  As  previously  mentioned  (see 
Section  3.2.2),  the  envelope  detected  output  of  the  bandpass  filter  has  been  shown  by  Rice 
[14]  to  have  a  Rayleigh  PDF.  Therefore,  the  PDF  for  coherent  detection  with  a  speckle 
target,  p,.j,  is  given  as  [10] 


Pcs(0  ~  ■ 


2/ 


{^CoS.SN  )  {k'oh.SN  )  ■*"  {^Diffusf  ) 


(3.3.12) 


where  i  represent  the  instantaneous  peak  envelope  detected  signal..  This  same  distribution 
could  be  found  by  substituting  zero  in  for  in  Eq.  (3.3.11).  An  example  of  the 
Rayleigh  PDF  is  shown  in  Fig.  6  where  the  mean  squared  values  used  are  shown. 
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Figure  6:  This  shows  a  Rayleigh  probability  density  with  an  increasing  diffuse 
component. 


3.3.22  A  Glint! Specular  Target  produces  a  deterministic  (non-statistically  varying)  return 
[2].  The  combination  of  this  return  signal  and  the  Gaussian  noise  gives  a  complex 
Gaussian  distribution  for  the  overall  detected  current.  The  return  is  complex  because 
there  are  both  in-phase  and  out-of-phase  components  due  to  differences  in  range  to  the 
target.  Tiie  PDF  representing  the  envelope  of  the  current  fluctuations  (signal  and  noise)  at 
the  output  of  the  bandpass  filter  for  coherent  detection  with  a  glint  target,  p,.f. ,  is  given  by 


[1,2,12,14] 


(3.3.13) 


This  is  known  as  the  Rician  PDF.  This  same  distribution  can  be  obtained  by  substitution 
in  zero  for  {jl^)  in  Eq.  (3.3.1 1)  since  this  discussion  is  for  an  entirely  glint  target.  An 

example  of  different  Rician  distributions  with  no  diffuse  component,  a  mean  glint 
component  of  10  nA,  and  an  increasing  noise  component  are  shown  in  Fig.  7. 
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3.4  The  Probabilities  of  False  Alarm 

A  false  alarm  occurs  when  a  return  signal  has  been  declared  a  target  when  in 
actuality  there  is  no  target  present.  For  a  given  threshold  level,  the  probability  of  false 
alarm  is  the  probability  that  the  noise  level  will  exceed  that  threshold  level. 
Mathematically,  it  is  defined  as  the  area  under  the  noise  PDF  which  exceeds  the  set 
threshold  level,  ij.  To  find  fhe  probability  of  false  alarm,  the  area  under  the  noise 
probability  density  curve  is  calculated  from  the  threshold  level  to  infinity  (see  Figure  1). 
As  seen  in  the  previous  sections,  the  dominating  noise  for  each  detection  technique  is 
different,  therefore  there  are  two  different  probability  density  functions  representing  the 
noise.  The  following  two  sections  discuss  the  probability  of  false  alarm  for  incoherent  and 
coherent  detection. 
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3.4. 1  The  Probability  of  False  Alarm  for  Incoherent  detection 

As  stated  above  the  probability  of  false  alarm  is  the  area  under  the  probability 
density  function  of  the  noise,  greater  than  a  set  threshold  level.  For  incoherent  detection 
the  probability  distribution  for  the  signal  shot  noise  current  fluctuation  is  a  Gaussian  PDF 
given  in  Eq.  (3.2.1)  as  [2,12] 


(3.4.1) 


where  /  is  the  instantaneous  detector  output  current  and  mean  squared 

noise  current  for  incoherent  detection.  The  probability  of  false  alarm  can  then  be 
expressed  as 


where  ij  is  the  threshold  current  level.  Equation  (3.4.2)  can  be  simplified  for  easier  use  by 
substituting  in  the  complementary  error  function,  ERFC(x),  where 


ERFC{x)  =  -^je-''dt. 


(3.4.3) 


Substituting  Eq.(3.4.3)  into  Eq.  (3.4.2)  results  in 


(3.4.4) 


3.4.2  The  Probability  of  False  Alarm  for  Coherent  Detection 


The  noise  probability  density  funcdon  for  coherent  detection  was  given  in  Eq. 
(3.2.2)  as  [12] 


Pco{i)  =  -r— 

ycohi 

The  probability  of  false  alarm  is  therefore 

Pcop^  =  J  Pco{i)di 


(  .-2  'I 

1 

rexp 

/ 

Pcop^  =  exp 


i 

1 

J  77. 

2  \  “P 

if  ycoh,sN  /  1 

i  ' 

#»vn! 

f  -JL 

J 

1 

\di 


27/2  T 

^yCoh,SN  / 


(3.4.5) 


(3.4.6) 


A  desired  probability  of  false  alarm  can  be  obtained  by  knowing  the  mean  squared 
noise  current  and  by  then  manipulating  the  threshold  current  if.  Once  the  desired 
probability  of  false  alarm  is  achieved  the  threshold  current  used  to  achieve  that  probability 
of  false  alarm  can  be  used  to  determine  the  probability  of  detection. 


3.5  The  Probability  of  Detection 

The  probability  of  detection  is  the  probability  that  a  received  signal  is  correctly 
declared  a  return  from  a  target.  Mathematically  speaking,  the  probability  of  detection  is 
the  area  under  the  probability  density  function  representing  the  signal  and  noise  greater 
than  the  decision  threshold  (see  Figure  1).  Since  there  are  two  detection  techniques  of 
interest  and  two  targets  of  inteiest  the  following  discussion  is  broken  down  into  four 
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sections.  Each  section  presents  the  probability  of  detection  for  a  specific  detection  scheme 
for  an  individual  type  of  target. 

3.5.1  The  Probability  of  Incoherent  Detection 

The  probability  of  detection  for  incoherent  detection  will  be  discussed  in  this 
section.  To  find  the  probability  of  detection  the  probability  density  function  for  each  type 
of  target  will  be  integrated  from  a  threshold  level,  determined  by  calculating  the 
probability  of  false  alarm,  to  infinity. 


3. 5. 1.1  The  Probability  of  Incoherent  Detection  with  a  Glint  Target  is  the  integral  of  the 
probability  density  function  for  combined  signal  and  noise.  This  probability  density 
function  was  given  in  Eq.  (3.3.4)  as  [12] 


P/c(0  = 


1 

(i-{iaunt)y 

P^iiiLoKSN  )  + 

^{{hncoh,SN  )  )) 

(3.5.1) 


The  probability  of  detection,  PdjQ,  is  then  calculated  as 


PdiG  — 


hncoKSN  )  )) 


Jexp 


^{{hncokSN  )  )) 


Vli.  (3.5.2) 


This  can  be  simplified  by  expressing  the  equation  in  terms  of  the  complementary  error 
function  (ERFC),  when  (<V  function  (ERF)  when 

('V  -  )  <  O) .  We  thus  find 
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O.SERFC  ^  where  i, -{i„„)  £  0 

^  ( ^{{hKoh,SN  )  ■*■  )))  ^ 

r  r  /  \2  V 

0.5  H-£«f  i{‘o«)-h) -  where <0 


(3.5.3) 


where  the  error  function  is  represented  as 


ERF{z)=-^\e-‘'dt, 

V^’  0 


(3.5.4) 


and  if  is  the  threshold  c.urent  discussed  earlier.  Once  ij  is  determined  using  the 
probability  of  false  alarm  the  probability  of  detection  can  be  determined. 


3.5.12  The  Probability  of  Incoherent  Detection  for  a  Speckle  Target  is  the  integral  from 
a  threshold  to  infinity  over  the  probability  density  function  given  in  Eq.  (3.3.10)  which  is 


{k-(,{n)  +  hW)f 


exp(-fliy) 

r(M) 


(3.5.5) 


So  the  probability  of  detection  is 


{k-{{n)+bW)f]a^W^-^  exp(-alV) 

2{{n)-^bW)  J  r(M) 


(3.5.6) 


where  xj  is  the  threshold  current  expressed  in  terms  of  photoelectrons  (Photoelectrons 
can  be  converted  to  current  by  multiplying  by  q//it  where  q  is  the  charge  on  an  electron 
and  At  is  the  transmitted  pulse  length  which  will  be  discussed  later  in  Section  4.1.2.)  This 
is  somewhat  of  an  intimidating  expression.  To  simplify  this  expression  the  outer  integral 
can  be  reduced  using  the  ERFC(z)  function  given  in  Eq.  (3.5.4),  ITie  result  is 
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m 

Fd,s=ijE/tFC 
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Xt  -‘((n)+b]V))  exp(-atV) 

^2((/i>W)  J  ■  r(Af) 


(3.5.7) 


3.5.2  Probability  of  Coherent  Detection 

In  the  following  section  the  probability  of  detection  will  be  discussed  for  coherent 
detection  with  a  glint  and  a  speckle  target. 


3.52.1  The  Probability  of  Coherent  Detection  with  a  Glint  Target  is  calculated  using  the 
Rician  PDF  given  in  Eq.  (3.3.13)  as 


(3.5.8) 


The  probability  of  detection  for  coherent  detection  with  a  glint  target,  Pd^Q:  is  then  given 

by 


where  due  to  the  Bessel  function  in  the  equation  there  is  no  clear  way  to  simplify  this 
equation  further. 


3. 5. 2. 2  The  Probability  of  Coherent  Detection  with  a  Speckle  Target  is  the  integral  from 
a  threshold  current  to  infuiity  over  a  Rayleigh  probability  distribution  given  in  Eq. 
(3.3. 12).  The  probability  of  detection  for  this  case  is 
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which  can  be  easily  integrated  and  reduced  to 


(3.5.10) 


(3.5.11) 


Now  that  the  equations  for  the  probability  of  detection  are  known  the 
characteristics  of  the  experimental  setup  need  to  be  known  in  order  to  do  a  numerical 
analysis  to  determine  the  probability  of  detection  for  each  case  for  both  detection  schemes. 


CHAPTER  IV 


4.0  SYSTEM  CONFir.llRATlON 


This  chapter  consists  of  a  detailed  description  of  the  components  of  the  LADAR 
system  used  to  compare  coherent  to  incoherent  detection.  This  discussion  consists  of  a 
component  by  component  description  of  the  system  and  how  each  component  aids  in  the 
operation  of  the  system  and  a  discussion  on  the  optimization  of  the  both  detection 
schemes. 


4.1  System  Layout 

4.1.1  Master  Oscillator 

A  top  view  of  the  laser  radar  system  is  shown  in  Figure  8.  The  master  oscillator 
(MO)  in  the  system  is  a  diode  pumped,  CW,  room  temperature,  Thulium  (Tm),  Holmium 
(Ho):YAG  laser.  This  laser  was  developed  for  this  project  by  CLR  Photonics  Inc.  This 
laser  has  a  maximum  output  of  80  mW.  The  output  of  the  MO  is  used  for  injection 
seeding  the  slave  oscillator  and  it  is  also  used  as  the  local  oscillator  (LO)  signal  when  the 
ladar  system  is  used  in  the  coherent  mode  of  detection.  The  first  mirror  and  most  of  the 
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Figure  8:  Layout  of  LADAR  system. 


turning  mirrors  in  the  system  are  coated  to  be  99.9%  reflective  at  45°.  The  next  element 
in  the  system  is  a  Faraday  optical  isolator,  which  was  acquired  from  Optics  for  Research, 
Model  #IO-4-HoYAG.  This  element  consists  of  an  input  polarizer,  a  rotator  material  and 
an  analyzer.  The  input  polarizer  insures  that  the  incoming  light  is  linearly  polarized 
horizontally.  The  rotator  material  then  rotates  the  polarization  of  the  light  by  45°  and  the 
beam  exits  through  the  analyzer,  which  is  positioned  at  45°.  Any  back-reflections  from 
other  components  in  the  system  enter  through  the  analyzer  and  the  polarization  is  rotated 
by  the  rotater  material  and  the  polarization  will  be  perpendicular  to  the  final  polarizer. 
This  gives  an  extinction  of  approximately  40dB.  The  optical  transmission  efficiency  of  this 
device  is  50%. 

The  next  element  in  the  system  is  a  half-wave  plate  which  is  used  to  control  the 
polarization  of  the  light  that  enters  the  slave  oscillator  (SO).  The  following  element  in  the 
system  is  a  98%  beam  splitter.  The  reflected  98%  is  used  as  the  LO  when  the  system  Ls 
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being  used  in  the  coherent  detection  scheme.  This  beam  is  reflected  around  the  end  of  the 
SO  and  reflected  so  the  beam  is  traveling  parallel  behind  the  cavity.  The  beam  then 
encounters  a  60mm  and  120mm  lens  pair  that  expands  the  beam  from  2mm  to  4mm.  The 
combination  of  the  next  two  elements  in  the  system,  the  half  wave  plate  and  the  thin  film 
Brewster's  angle  polarizer,  allow  control  of  the  amount  of  light  incident  on  the  detector 
which  is  used  as  the  LO.  The  thin  film  Brewster's  angle  polarizer  was  obtained  from 
Rocky  Mountain  Instruments.  The  Brewster's  angle  for  2.09|xm  for  this  element  is  57°. 
The  light  that  reflects  off  the  thin  film  Brewster's  angle  polarizer  is  then  reflected  into  the 
detection  optics  by  a  15%  beam  splitter.  The  detection  optics  will  be  discussed  later. 

The  transmitted  2%  from  the  98%  beam  splitter  is  frequency  shifted  27. 1  MHz, 
into  the  first  order,  by  an  acoustic  optic  modulator  (AOM).  This  AOM  was  purchased 
from  Newport  Electro-Optical  Systems,  Model  #N31027-4D.  The  diffraction  efficiency 
of  this  AOM  for  2.09p.m  is  about  60%.  The  beam  that  is  frequency  shifted  into  the  first 
order  is  then  reflected  into  the  SO  for  injection  seeding. 

4.1,2  Slave  Oscillator 

The  slave  oscillator  is  a  Chromium  (Cr),  Tm,  Ho:YAG,  flaslilamp  pumped,  Q- 
switched  laser.  The  end  mirrors  are  greater  than  99.5%  reflective  at  0°  incidence,  planar, 
AR  coated  at  2,09pm  and  have  a  damage  threshold  of  lOOMW/cm^,  The  output  coupling 
is  controlled  by  the  combination  of  a  quarter-wave  plate,  the  end  mirror  and  a  thin  film 
Brewster's  angle  polarizer. 

The  polarization  of  the  light  in  the  cavity  is  horizontal  until  it  transmits  through  the 
quarter  wave  plate  as  seen  in  Fig,  9.  Once  the  light  translates  through  the  quarter-wave 
plate  the  polarization  of  the  tight  is  generally  elliptical.  ITie  rotation  of  this  ellipitically 
polarized  light  reverses  direction  after  tf>e  light  strikes  the  end  mirror.  Once  the  light 
passes  through  the  quarter-wave  plate,  on  its  return  trip,  the  light  will  again  be  linearly 
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polarized  but,  there  will  be  both  a  horizontal  and  a  vertical  component.  The  vertical 
component  of  the  polarization  is  reflected  out  of  the  cavity  when  it  strikes  the  thin  film 
Brewster's  angle  polarizer.  The  amount  of  vertically  polarized  light  depends  on  the 
position  of  the  fast  axis  of  the  quarter-wave  plate.  The  quarter-wave  plate  is  set  so  the 
output  coupling  of  the  cavity  is  20%.  When  the  output  coupling  is  set  lower  the  intensity 
of  the  light  in  the  cavity  is  such  that  optical  elements  can  be  and  have  been  damaged. 

The  Kigre  Inc.  FOM-455C  laser  pump  cavity  is  dual-lamp  with  a  diffusely- 
reflecting  elliptical  cavity.  The  diffusely-reflecting  cavity  is  less  efficient  than  a  reflector 
cavity,  but  the  pumping  is  more  uniform,  giving  a  more  uniform  output  beam  and  better 
quality  thermal  lensing  characteristics.  Thermal  leasing  occurs  in  the  laser  rod  because  of 
heating  from  the  flashlamps  and  the  cooling  from  the  water.  The  rod  is  water  cooled  but 
the  rate  at  which  the  rod  is  pumped  determines  how  cool  the  rod  is  kept.  As  the  repetition 
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rate  increases,  the  rod  begins  to  cool  unevenly  causing  a  temperature  gradient  between  the 
outside  and  inside  of  the  rod,  where  the  inside  of  the  rod  is  at  a  higher  temperature  than 
the  outside.  Tliis  temperature  gradient  causes  a  radial  change  in  the  index  of  refraction. 
This  change  in  the  index  causes  the  rod  to  act  like  a  lens,  hence  the  name  of  thermal 
lensing.  Wlien  thermal  lensing  occurs  the  mode  volume  decreases,  decreasing  the  output 
power  of  the  laser.  If  the  pump  rate  is  fast  enough  the  rod  does  not  have  enough  time  to 
cool  to  give  as  great  a  temperature  differential  and  it  will  only  decrease  the  mode  volume 
slightly  [17].  Since  the  mirrors  in  this  cavity  are  planar  the  cavity  is  critically  stable,  but 
with  thermal  lensing  occurring,  the  stability  of  the  cavity  increases  due  to  the  focusing 
power  of  the  rod. 

Thie  pump  chamber  is  a  barium-sulfate-packed  KC331  glass  filter/flow  tube.  The 
glass  filter  absorbs  energy  below  343  nm  which  are  not  needed  to  pump  the  crystal  to 
population  inversion.  This  glass  also  absorbs  heat  preventing  unnecessary  heating  of  the 
rod.  The  flow  tube  allows  the  Kigre  Xenon  flashlamps  and  the  rod  to  be  cooled  with 
deionized  water  to  20°C.  The  cooling  of  the  flashlamps  and  the  laser  rod  is  controlled  by 
a  Neslabs  HX  recirculating  chiller. 

Tlie  YAG  laser  rod  is  55  mm  long  and  4mm  in  diameter.  The  dopant 
concentrations  are  0.86%  Chromium  (Cr),  5.83%  Thulium  (Tm),  and  0.34%  Holmium 
(Ho).  The  energy  level  diagram  for  Cr,  Tm,  Ho:YAG  is  shown  in  Fig.  10.  The  flashlamp 
energy  is  absorbed  by  the  Cr^"*"  ion.  Once  the  atoms  .settle  in  the  ^T2  and  states  they 
transfer  to  the  ^^^4  ^  Tm3+  ion.  The  excited  Tm^"^  atom  then  interacts 

with  a  ground  state  Tm^"*"  in  a  cross-relaxation  process.  This  process  produces  two  Tm3+ 
ions  in  the  ^¥4.  state.  The  energy  is  then  transferred  to  the  ^17  state  where  lasing  occurs  in 
the  ^17  _  transition  [6]. 
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The  slave  oscillator  is  Q-switched  using  Fast  P'dse's  3902W-2100  LiNb03  Pockels 
cell.  The  Q  of  a  cavity  is  defined  as  the  energy  stored  in  a  system  at  resonance  divided  by 
the  energy  lost  in  a  cycle  of  oscilladon.  So  the  higher  the  Q  the  lower  the  losses.  Q- 


Figure  10:  Energy  level  diagram  for  Cr,Tm,Ho;  YAG,  laser  rod. 

switching  is  the  switching  between  a  cavity  of  a  low  Q,  which  prevents  lasing,  to  a  cavity 
with  a  high  Q,  which  allows  the  cavity  to  lase.  The  changing  of  the  Q  of  the  transmitter  Ls 
controlled  by  the  LiNb03  Pockels  cell  using  the  electro-optic  effect.  Enough  voltage  is 
applied  to  the  Pockels  cell  to  rotate  the  polarization  inside  the  cavity  so  that  the 
combination  of  the  quarter- wave  plate  and  the  thin  film  polarizer  reflects  the  light  out  of 
the  cavity  causing  it  not  to  lase.  T>-":  voltage  that  would  normally  be  applied  to  the 
Pockels  cell  to  accomplish  this  task  is  the  quarter-wave  voltage.  The  quarter-wave 
voltage  for  this  cell  is  2.5  kV,  but  the  cavity  of  this  laser  is  veiy  lossy  so  a  voltage  much 
lower  is  sufficient  to  perfomi  the  Q-switching.  Experimentally  this  voltage  has  been  found 
to  be  about  960  Volts.  It  was  also  found  that  this  voltage  is  directly  dependent  on  both 
the  repetition  rate  and  the  pump  energy.  An  example  of  an  unQ-switched  puise  is  shown 
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in  Fig.  11  and  a  Q-switched  pulse  is  shown  in  Fig.  12.  The  two  main  differences  of 
concern  are  the  length  of  the  pulse  and  the  magnitude  of  the  pulse.  The  unQ-switched 
pulse  is  about  2  ms  in  length  whereas  the  Q~switched  pulse  is  about  1  ms  in  length  with  a 
FWHM  of  500ns.  The  vertical  scales  in  both  Figs.  1 1  and  12  are  in  volts,  which  was  an 
observed  output  from  the  InGaAs  detector.  The  output  of  the  Q-switched  pulse  is  nearly 
an  order  of  magnitude  greater. 


Figure  11:  Slave  oscillator  long  pulse  output. 

Another  element  in  the  SO  cavity  is  a  75|J.m  etalon.  This  etalon  is  used  to  tune  the 
wavelength  of  the  cavity.  By  changing  the  angle  of  the  etalon  in  the  cavity  the  wavelength 
can  be  changed. 

In  order  to  do  coherent  detection  the  wavelength  of  the  SO  is  locked  to  that  of  the 
MO  except  for  the  27.1  MHz  frequency  shift  introduced  by  the  AOM.  In  order  to 
properly  injection  seed,  the  length  of  the  SO  cavity  has  to  be  such  that  the  desired 
wavelength  will  oscillate  in  the  cavity.  To  control  this  length,  one  end  mirror  is  mounted 
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to  a  piezoelectric  transducer  allowing  it  to  be  dithered  to  the  correct  wavelength  to  lase. 
A  servo- loop  was  created  by  CTI  in  Boulder  Colorado  to  do  this  job.  There  is  a  InGaAs 
photo-conductor,  purchased  from  Epitaxx,  mounted  in  the  end  of  the  laser  cavity.  This 


Figure  12:  Slave  oscillator  output  when  Q-switched. 

detector  sends  its  signal  to  a  servo  box  which  monitors  the  build-up  time  of  the  la.ser 
pulse.  When  the  laser  is  properly  injection- seeded  the  laser  pulse  will  occur  earlier  than 
normal.  What  this  servo  loop  does  is  monitor  the  build-up  time  of  the  pulse  and  it  dithers 
the  cavity  to  decrease  this  build-up  time.  There  is  a  threshold  level  that  can  be  tuned  on 
the  servo  system  so  an  alarm  will  go  off  when  the  system  is  not  injection- seeding. 

Even  though  there  is  a  servo-loop  to  aid  in  injection-seeding  the  SO  and  the  MO 
have  to  be  lasing  at  nearly  the  same  wavelength  and  the  MO  has  to  be  aligned  so  that  it 
oscillates  in  the  cavity  along  the  same  path  that  the  SO  lasts.  The  etalon  in  the  .slave 
cavity  allows  the  SO  to  be  wavelength  tuned  as  does  the  etalons  in  the  MO.  The 
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wavelength  that  the  system  is  tuned  to  is  2.09046|im.  Both  the  MO  and  the  SO  have  to 
be  close  to  this  wavelength  or  the  system  will  not  injection  seed. 

The  SO  laser  requires  900V  applied  to  the  flashlamps,  which  results  in  about  96  J 
of  applied  energy.  The  repetition  rate  is  2.7  Hz  with  a  long  pulse  output  of  20  mJ.  The 
shape  of  the  long  pulse  is  seen  if  Fig.  1 1.  When  the  LiNb03  Q-switch  is  used  the  output 
reduces  to  approximately  17  inJ  with  a  FWHM  pulse  length  of  500  ns.  An  example  of  this 
pulse  is  seen  in  Fig.  12. 

The  output  pulse  fiom  the  transmitter  that  reflects  out  of  the  cavity,  is  first 
expanded  2x  by  a  50mm  and  a  lOOmm  lens  pair.  The  beam  then  reflects  off  of  two 
mirrors  that  are  used  to  walk  the  beam  during  alignment  procedures.  The  beam  then 
strikes  another  thin  film  Brewster's  angle  polarizer.  Since  this  beam  is  vertically  polarized, 
it  reflects  off  of  the  polarizer.  There  is  some  leakage  of  the  outgoing  pulse  through  the 
polarizer.  This  portion  is  used  to  monitor  the  outgoing  pulse  to  the  target  using  a  JDIOOO 
Molcctron  energy  meter.  The  ratio  of  the  light  reflected  to  that  which  is  transmitted 
through  the  polarizer  has  been  measured  as  15.3  to  1. 

The  thin  film  Brewster's  angle  polarizer  and  the  quarter-wave  plate  make  up  the 
transmit  receive  switch  for  the  laser  radar  system.  The  beam  reflects  off  the  polarizer 
transmits  through  the  quarter  wave  plate,  which  produces  circular  polarization.  The 
direction  of  the  circular  polarization  is  reversed  when  the  beam  strikes  a  target.  So  when 
the  return  beam  travels  through  the  quarter  wave  plate  the  light  is  changed  from  circular 
polarization  to  linear  (horizontal)  polarization.  Being  horizontally  polarized  the  light 
transmits  through  the  polarizer  and  the  polarization  is  then  rotated,  using  a  half  wave 
plate,  to  match  that  of  the  LO  beam  incident  on  the  detector. 

Before  the  beam  is  transmitted  to  the  target  it  is  expanded  20x  by  a  Lambda/lO's 
Laser  Beam  Expander  Model  LBX.  This  beam  expander  is  an  afocal,  decentered,  Dall- 
Kirkhan  telescope.  The  primary  mirror  is  a  decentered  section  of  a  concave  ellipsoid  while 
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the  secondary  mirror  is  a  convex  sphere.  The  beam  leaving  the  telescope  is  steered  using 
a  mechanical  beam  steerer.  The  beam  steering  device  has  a  6  inch  flat  mirror  mounted  on 
a  two  axis,  computer  controlled,  Arotech  gimble. 


4.2  Detection  Scheme 

The  theoretical  an.Uysis  has  been  done  assuming  optimum  detector  performance 
for  each  detection  scheme.  For  incoherent  detection  the  optimum  detector  performance 
occurs  with  the  best  filling  of  the  detector  with  the  incoming  light.  For  coherent  detection 
the  optimum  detection  of  the  light  is  dependent  on  the  efficiency  of  the  mixing  between 
the  local  oscillator  and  the  return  pulse.  Since  this  analysis  has  not  taken  into  account 
turbulence  the  data  taken  needs  to  be  taken  together  or  one  set  directly  after  the  other. 
This  allows  the  data  to  be  taken  with  nearly  the  same  atmospheric  conditions.  To  save 
time  transitioning  betv/een  detection  schemes,  a  mirror  on  a  translation  stage  has  been 
placed  in  the  return  beam  path.  With  the  translation  of  the  mirror  into  or  out  of  the 
receiver  path,  the  return  energy  can  be  sent  to  either  detection  scheme . 

For  incoherent  detection  the  return  light  that  transmits  through  the  15%  beam 
splitter,  a  38  nm  optical  bandpass  filter  from  Optical  Filter  Corporation  and  is  focused 
down  onto  an  InGaAs  detector.  The  return  energy  from  the  target  is  returning  as  a  plane 
wave.  Since  it  is  a  plane  wave  the  focused  return  energy  onto  the  detector  is  represented 
as  an  Airy  disk  [18].  The  intensity  pattern  of  the  Airy  disk  is  shown  in  Fig.  13,  which  is 
mathematically  represented  as 


(4.2.1) 
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where  J](x)  is  a  first  order  Bessel  function,  and  x^k^aq^jR^,  where  a  is  the  limiting 
optical  aperture,  q^  is  the  radius  of  tlie  airy  disk,  is  the  wave  number  and  Ri  is  the 
distance  from  the  aperture  to  the  image.  For  this  system  it  was  decided  that  for  incoherent 
detection  both  the  central  lobe  and  the  first  ring  will  be  used.  The  first  ring  is  zero  at 
j:=7.016.  So  using  this  and  the  equation  for  jc  the  value  for  the  focal  length  was  found  to 
be  53mm.  The  50  mm  focal  lens  was  already  being  used  to  expand  one  of  the  beams  in 
the  system  so  a  40  mm  focal  length  lens  was  substituted.. 


Figure  13:  One  side  of  Airy  disk  pattern. 


For  coherent  detection  the  optimum  detection  of  the  light  occurs  when  only  the 
first  lobe  of  the  Airy  disk  pattern  for  the  combination  of  the  local  oscillator  and  the  return 
signal  is  incident  bn  the  detector.  The  first  ring  is  180°  out  of  phase  from  the  centra]  lobe, 
so  mixing  the  first  ring  with  the  central  lobe  would  decrease  tlie  heterodyne  efficiency. 

The  heterodyne  efficiency  is  a  ratio  that  expresses  how  efficient  the  returii  signal 
and  the  LO  mix  spatially.  When  the  return  signal  and  the  LO  are  matched  Airy  functions 
at  the  plane  of  the  detector,  the  heterodyne  efficiency,  7  ,  is  given  by  [18, 19] 


(4.2.2) 
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where  y  =  MrD^j  fk^fis  the  focal  length  of  lens,  £)g  is  the  diameter  of  the  laser  beam,  X,  is 
the  wavelength,  r  is  the  radius  of  the  detector,  and  and  J^  are  Bessel  functions  of  the 
first  kind.  Equation  4.2.2  is  plotted  in  Fig.  14. 

For  our  system,  the  radius  of  the  detector  r  is  50  pm,  X  is  2.09  pm,  and  the 
diameter  of  the  laser  beam,  Dg,  is  4  mm.  Notice  now  thi't  the  smaller  the  focal  length  lens 


Figure  14:  Theoretical  heterodyne  efficiency. 

chosen,  the  better  the  heterodyne  efficiency.  For  our  work,  a  reasonable  focal  length  lens 
of  80  rnm  was  chosen,  giving  an  y  of  3.75  and  a  theoretical  heterodyne  efficiency  of  0.84. 

A  schematic  of  the  detection  packages  provided  by  Coherent  Technologies,  Inc. 
(CTI)  of  Boulder,  CO  is  shown  in  Fig.  15  [20].  The  photudetector  is  reversed  biased, 
Vp,  by  a  3  Volt  battery  and  is  in  series  with  a  1  kO  resistor.  The  photodetector  used  is  a 
p-i-n  InGaAs  photodiode  which  has  a  measured  daik  cunent  of  135  nA  and  responsivity 
of  l.l  A/W  at  2.09  pm.  Since  for  coherent  detection  the  detected  signal  will  be  at  high 
frequencies,  27.1  MHz,  and  for  direct  detection  the  detected  signal  will  be  at  low 
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frequencies,  2.7  Hz,  the  post  detection  electronics  are  different  for  each  detection  scheme. 
For  incoherent  detection  the  amplifier  was  purchased  from  Analog  Modules  and  it  has  a 
bandwidth  of  35  MHz  and  a  midband  gain  of  60  dB  [21],  There  are  no  filters  used  for  the 
incoherent  detection  scheme  so  the  bandwidth  is  limited  by  the  amplifier.  The  optimum 
load  resistance,  Rl,  found  by  CTI  for  incoherent  detection  is  16  Mi2  [20]. 
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Figure  IS:  Detection  package  from  CTI. 

For  the  heterodyne  detection  case,  the  same  reverse  bias  voltage,  series  resistance 
and  photodiode  are  used  (see  Figure  15).  The  amplifier  is  different,  being  provided  by 
Miteq,  it  has  a  bandwidth  of  100  MHz  and  a  midband  gain  of  53  dB  [22].  To  limit  the 
bandwidth  of  the  detection  package  and  to  isolate  the  return  signal,  a  41  MHz  bandpass 
filter  was  used  centered  at  the  intermediate  fiequency,  27.1  MHz  [20].  The  decrease  in 
the  signal  due  to  the  bandpass  filter  is  less  than  1  dB.  The  load  resistance  found  by 
Cn  for  optimized  heterodyne  detection  for  this  system  at  27. 1  MHz  is  450D  [20]. 


CHAPTER  V 


5.0  THEORETICAI,  COMPUTATIONS 

This  chapter  will  use  the  theory  discussed  in  CHAPl'ER  III  and  the  values 
specified  for  the  components  in  the  detection  electronics  to  calculate  theoretical 
probabilities  of  detection  for  each  detection  scheme.  In  order  to  find  the  probability  of 
detection  the  limiting  mean  squared  noise  current  for  each  detection  scheme  needs  to  be 
identified.  The  limiting  mean  squared  noise  current  will  consist  of  a  sum  of  all  the  noi.se 
currents.  Some  of  the  these  noise  currents  will  be  small  enough  that  their  contribution  is 
miniiTial  and  can  be  ignored.  Once  the  sum  of  the  dominating  noise  currents  is  found  it 
will  be  used  in  the  noise  PDF  to  find  the  thresholds  needed  for  the  desired  probabilities  of 
false  alarm.  Once  the  thresholds  are  found,  they  can  be  used  along  with  the  limiting  mean 
squared  noise  current  in  the  probability  density  function  of  the  current  fluctuation  to 
calculate  the  probability  of  detection.  The  probability  of  detection  will  then  be  plotted 
against  signal-to-noise  ratio  to  produce  curves  for  each  probability  of  false  alarm. 

5.1  Evaluation  of  Noise  in  the  2  am  LADAR  System 

The  noise  sources  discussed  in  CHAPTER  III  are  dark  current,  Johnson, 
background,  shot,  and  amplifier  noise.  The  average  theoretical  values  for  these  noise 
sources  will  be  calculated  for  each  detection  scheme  using  the  known  characteristics  of  the 
system.  Once  the  magnitude  of  these  sources  is  know,  a  threshold  can  be  calculated 
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for  given  probabilities  of  false  alarm,  and  this  threshold  can  be  used  to  calculate  the 
probability  of  detection  for  each  case. 


S.1.1  Dark  Current  Noise 

Dark  current  shot  noise  is  directly  proportional  to  the  dark  current  for  a  particular 
detector  and  each  detector  is  unique  even  if  it  comes  from  the  same  processing  lot.  The 
dark  current  for  one  of  the  purchased  p-i-n  strained  InGaAs  photodiodes  was  measured  to 
be  13S  nA  for  the  3V  reverse  bias  used  in  the  detection  bias  circuit.  For  the  coherent 
detection  scheme  the  bandwidth  is  determined  by  the  bandpass  filter  to  be  41  MHz  [22] 
This  gives  a  mean-square  dark  current  shot  noise  for  coherent  detection  of 

W*  ~^^^Dark^ 

=  2(1.602  X  10-'’C)(135  x  10"’  A)(41  x  10*  Hz)  (5.1.1) 

=  1.77x10"'* 

The  only  difference  between  the  two  detection  schemes  is  the  bandwidth  and  for 
the  incoherent  detection  scheme  the  Analog  Modules  amplifier  determines  the  bandwidth. 
The  bandwidth  of  the  amplifier  is  35  MHz  [21].  This  gives  a  mean  square  dark  current 
shot  noise  for  incoherent  detection  of  1.51x10'^^  A^. 


5.1.2  Johnson  Noise 

Johnson  noise  is  inversely  proportional  to  the  load  resistance.  According  to  CTTI 
the  optimum  load  resistance  for  coherent  detection  is  450f2  [20].  Therefore,  the  Johnson 
noise  for  the  coherent  case  is 


4^rg 

^Load 

4(l.38  X 10"“  /^)(293A')(41  X 10* 
450ft 


=  1.47xlO-‘*A^ 


(5.1.2) 
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For  incoherent  detection  CTl  gave  an  optimum  load  resistance  of  The 

mean  squared  Johnson  noise  current  for  the  incoherent  case  is  then  3.54x1 This  Ls 
much  smaller  than  that  for  the  coherent  case  but  that  was  expected  with  the  larger  load 
resistor. 

5.1.3  Background  Noise 

For  this  analysis  the  worst  case  scenario  for  background  noise  is  assumed,  which 
is,  the  scenario  of  the  LADAR  looking  at  a  sun  illuminated  cloud.  The  amount  of  solar 
power  received  at  the  detector  was  given  in  CHAPTER  III  as 

^SB  (5.1.3) 

where  it,  is  the  fraction  of  solar  radiation  penetrating  the  Earth’s  atmosphere,  Sjrr  is  the 
solar  irradiance,  AX  is  the  wavelength  band  of  the  optical  band  pass  filter  centered  at  2.09 
^m  (see  chapter  4),  is  the  solid  angle  over  which  the  receiver  collects  energy  from  the 
radiating  body  (easily  found  by  dividing  the  area  of  the  beam  at  the  target  by  the  range  to 
the  target  squared),  p  is  the  target  reflectivity,  rjsys  is  the  ladar  system  optical  efficiency 
which  includes  the  reflectivity  of  the  mirrors  in  the  system  and  all  beam  splitters,  and  A^  is 
the  receiver  area.  The  values  used  for  these  variables  and  their  source  are  included  in 
Table  2. 

Substituting  the  values  in  Table  2  into  (5.1.3) 

=  (0. 95)(o.  0 1 380  A  J(2. 1 25  X 1 0^  jr  )(1)(0. 65)(0. 00785) 

=  3.91xl0"‘V 


(5.1.4) 
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Table  2: 


Summary  of  values  for  backgrousid  noise  calculation. 


Variable 

Value 

Units 

Source  of  Infomnation 

^’i 

0.95 

RCA  Electro-Optics  Handbook,  [23] 

SiRR 

0.01 

W/m2A 

It 

AX 

380 

A 

Optical  Filter  Corporation(OFC)  [24] 

Q/f 

2.125x1 

steradians 

Area  of  Propagated  Gaussian  beain/Range2,  range=20()0m 

P 

0’ 

Worst  case  scenario. 

^SYS 

1 

Actual  measurement  of  the  optical  system  tluxiughput 

Ar 

0.65 

.00785 

m2 

Receiver  area  {)f  10  cm  telescop:^ 

The  shot  noise  current  produced  at  the  detector  by  the  background  power  is  [4], 

=  (5.1.5) 

where  Rp„  is  the  detector  responsivity  in  Amperes/Watt,  q  is  the  charge  on  an  electron  in 
Coulombs,  and  B  is  the  electrical  bandwidth  of  the  receiver  in  Hertz.  Substituting  in  the 
appropriate  parameters,  we  find, 

(i)  =  2(l.602xl0-'’cX41xir/«2)^i.|A^3.9l5xl0-"lv)  ^ 

=  5.66xl0“'"/4\ 

For  incoherent  detection  the  bandwidth  is  35  MHz  so  the  background  noise  is  4.83x1 
a2. 


5.1.4  Amplifier  Noise 

Amplifier  noise  is  noise  added  to  the  signal  when  it  goes  through  the  amplifier. 
Amplifier  noise  is  a  thermal  noise  source  and  typically  the  noise  for  a  particular  amplifier  is 
specified  by  the  manufacture  using  an  equivalent  temperature,  Tj\  The  amplifier  obtained 
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for  this  systeirt  is  a  53  dB  amplifier  from  Mitcq,  model  AU-3A-01 10  and  the  typical  noise 
temperature  is  given  by  the  manufacturer  to  be  78  Kelvin.  The  amplifiers,  both  for  the 
coherent  and  the  incoherent  detection,  see  50il  for  the  load  resistance.  For  coherent 
detection  the  amplifier  noise  is 


R 


4(1. 38  X 10"“  yJ(78K^)(41  x  lO'  Hz) 


50fi 


(/Lp)  =  3.53x10-''A^ 


(5.1.7) 


For  the  incoherent  ca.se  an  equivalent  noise  temperature  was  not  given  by  the 
manufactuier  but,  a  measured  input  spectral  noise  density  was  given.  The  input  spectral 
noise  density  is  a  measure  of  the  noise  added  by  the  amplifier  but  it  is  specified  at  the  input 
to  the  amplifier.  Since  all  of  our  noise  sources  are  being  considered  at  the  input  of  the 

ainplifier  this  is  convenient.  The  spectral  noise  density  is  specified  by  the  manufacture  to 
be  0.612^.  In  order  to  get  the  noise  voltage  squared  al  the  input  of  the  amplifier  the 

spectral  noise  density  is  .squared  and  then  multiplied  by  the  bandwidth  as  seen  below. 

(vL,)  =  (0.612;^)'35x10"//z  =  1.31x10-"K"  (5.1.8) 


To  find  the  noise  current  before  the  amplifier  the  voltage  .squared  is  divided  by  the  input 
impedance  of  the  amplifier,  which  is  1  MQ. 


1.31x10~“K^ 

(ixio'n)^ 


1.31  X 10"“ 


(5.1.9) 


5./.5  Shot  Noise 

Shot  noise  is  (he  noise  produced  when  light  is  incident  on  a  detector.  This  noi.se  Ls 
present  in  both  coherent  and  incoherent  detection.  For  coherent  detection  the  desire  is  for 
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the  5»ystem  to  be  shot  noise  limited  by  illuminating  the  detection  with  the  local  oscillator 
beam.  For  the  incoherent  detection  scheme  shot  noise  is  produced  by  the  signal  return 
from  the  target.  In  this  section  shot  noise  will  be  calculated  for  both  cases. 


5. 1.5.1  Shot  Noise  in  the  Coherent  Detection  Case  is  produced  by  illuminating  the 
detector  with  the  local  oscillator.  To  insure  that  shot  noise  dominates  all  other  noise 
sources,  the  local  oscillator  power  incident  on  the  detector  was  increased  until  the  shot 
noise  was  10  dB  above  the  dominant  noise,  which  in  our  case  is  thermal  noise  (as  seen  in 
Section  5.1.2).  The  mean  squared  thermal  noise  current  was  calculated  to  be  1.47x10' 
a2  and  if  the  local  oscillator  is  increased  until  there  is  a  10  dB  difference  in  the  noise 
current  the  shot  noise  current  should  be  1.47xl0‘^4  A^.  These  noise  levels  can  be 
verified  for  the  system  on  the  spectrum  analyzer. 

Figure  16  shows  the  output  of  the  detection  electronics  displayed  on  the  spectrum 
analyzer  with  no  local  oscillator  incident  on  the  detector.  The  average  noise  level  seen  in 
Fig.  16  is  -1 17.5  dBm.  This  can  be  compared  to  the  dominating  noise  level  predicted  by 
finding  the  mean  squared  current  that  corresponds  to  this  noise  level.  To  do  this  the 
power  that  this  corresponds  to  can  be  calculated  as  shown  below. 


■117.5t/6m  =  101og 


f— 1 

\\mWj 


p  =  nnx\o-'^w 


(5.1.10) 


where  P  represents  the  amount  of  measured  power  with  a  10  Hz  resolution  bandwidth. 
To  correct  for  the  resolution  bandwidth  the  power  is  divided  by  the  10  Hz  and  then 
multiplied  by  the  41  MHz  bandwidth  of  the  detection  electronics.  This  gives  a  power 
measured  for  the  41  MHz  bandwidth  of  7.29x10'^  W.  This  was  measured  with  the  Miteq 
amplifier  in  place  which  has  a  53  dB  gain.  Since  all  the  noise  sources  have  been  calculated 
before  the  amplifier  the  effect  of  the  amplifier  needs  to  be  taken  into  account,  which  gives 
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Figure  16:  Spectrum  analyzer  display  of  output  of  coherent  detection  electronics  using  a 
10  Hz  resolution  bandwidth. 

a  power  of  3.65x10'^^  W.  Finally  we  want  the  mean  squared  current,  so  the  power  just 
calculated  needs  to  be  divided  by  the  load  resistance  to  get  the  mean  squared  current. 

>  =  ^  =  7.31X10-“  (5.1.11) 

The  dominating  noise  for  the  coherent  detection  case  is  the  sum  of  the  two  dominating 
noise  terms  thermal  and  amplifier  noise.  These  were  calculated  in  Eqs.  (5.1.2)  and  (5.1.7) 
and  the  sum  is  5.0x10'^^a2  which  is  a  factor  of  6.8  more  than  the  noise  measured 
experimentally.  A  factor  of  6.8  on  the  spectrum  analyzer  is  8  4  dB  and  the  spectrum 
analyzer  has  an  uncertainty  of  4  to  5  dB. 

The  master  oscillator  laser  was  then  used  to  illuminate  the  detector.  The  amount 
of  power  incident  on  the  detector  was  increased  until  there  was  a  10  dBm  difference  in  the 
amplitude,  which  was  1.25  mW  of  power.  This  result  can  be  seen  in  Fig.  17. 
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Figure  17;  Spectrum  analyzer  output  of  coherent  detection  electronics  with  1.25  mW  of 
power  incident  on  the  detector. 

The  average  noise  level  seen  in  Fig.  17  is  -108.3  dBm.  Using  the  same  analysis  as  used  for 
Fig.  16  the  mean  squared  shot  noise  current  is  6.08x10"^^  which  is  9.2  dB  less  than 
the  predicted  shot  noise  level  of  5.0x10" A^.  The  predicted  shot  noise  level  is  lOdB 
above  the  sum  of  the  thermal  and  amplifier  noise  terms.  The  lOdB  level  above  the  noise 
was  chosen  to  ensure  that  the  local  oscillator  would  dominate  the  noise. 

5. 1.5. 2  Shot  Noise  in  the  Incoherent  Case  is  produced  by  the  return  from  the  target  being 
incident  on  the  detector.  To  calculate  the  shot  noise  produced  in  the  incoherent  detection 
scheme  the  power  of  the  return  (Pr)  from  the  target  needs  to  be  determined,  once  the 
return  power  is  known  the  amount  of  shot  noise  can  be  calculated. 

To  find  the  return  power  the  radar  range  equation  is  used,  which  has  the  form  [4] 


iTyT  ^ 

r  niD  T 

— 

(5.1.12) 
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where  Pj  is  the  transmitted  power,  R  is  the  range,  Gg  is  the  beam  divergence,  Gj  is  the 
target  cross  section,  D  is  the  diameter  of  the  receiver,  and  t\atm  atmospheric 

transmission  factor.  The  target  cross  section  for  an  extended  target  is 

a7.=Jip/?^G8  (5.1.13) 

where  p  is  the  reflectivity  of  the  target, 
and  reduced  to 

The  values  used  for  these  variables  are  shown  in  Table  3. 

Table  3: 


For  the  transmitted  power  it  was  assumed  that  the  transmitted  pulse  was  10  mJ  and  the 
pulse  length  was  500  ns,  which  is  tj/pical  for  this  system.  The  range  chosen  is  typical  of 
those  used  during  the  experimental  process  and  target  reflectivity  was  chosen  earlier  to  be 
1  but  for  a  typical  lambertian  target  it  is  0.9.  The  atmospheric  attenuation  coefficient  was 
found  in  the  RCA  Handbook[25].  These  values  can  be  sub.stituted  into  Eq.  (5.1.14)  to 
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give  a  return  power  of  6.185  |xW.  This  can  be  use  in  the  equation  for  the  shot  noise  as 
shown  below. 

=  2(  1. 602  X  lO'*’  )i35MHz)(6. 1 85^1V)(  1. 1  (5.1.15) 

=  7. 63x10''’ A 

This  noise  happens  to  be  the  dominating  noise. 


Table  4: 


Noise  Source 

Detection  Scheme 

Coherent 

Incoherent 

Dark  current 

1.77x10- 1«a2 

1.51x10-18  a2 

Thermal 

1.47x10-15  a2 

3.54x10-20  a2 

Background 

5.66x10-22  a2 

4.83x10-22  a2 

Shot 

1.50x10-14  a2 

7.63x10-17a2 

Amplifier 

3.53x10-15  a2 

1.31x10-23  a2 

As  can  be  seen  from  Table  4,  coherent  detection  is  dominated  by  shot  noise  which 
is  greater  than  the  otherwise  dominating  noise  (the  sum  of  thermal  and  amplifier  noise) 
and  for  the  incoherent  detection  case  the  dominating  noise  is  the  shot  noise  produced  1^ 
the  return  signal  from  the  target.  This  seems  a  bit  surprising  but  the  target  used  was  only 
2  km  away,  and  it  is  not  quite  into  the  far  field  (far  field,  2.4  km).  In  order  to  find  the 
range  at  which  the  signal  shot  noise  no  longer  dominates,  Eqs.  (5.1.13)  through  (5.1.16) 
were  used  in  reverse  order  with  the  next  largest  noise  source,  which  is  dark  cunent  noise, 
to  find  the  range  at  which  shot  noise  from  the  signal  is  no  longer  the  doininating  noise. 
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This  back  calculation  was  done  and  the  range  at  which  signal  shot  noise  no  longer 
dominates  is  14.2  km.  This  is  probably  further  away  than  actually  needed  to  make  the 
sy.stem  non-signal  shot  noise  limited  because  the  atmospheric  transmission  decrea.ses  as 
the  range  increases. 


5.2  Calculation  of  the  Probability  of  False  Alarm 

The  major  noise  sources  were  analyzed  in  the  previous  section  and  the  values  of 
each  of  the  noi.se  .sources  was  determined.  U.sing  the  values  of  the  dominating  noi.se  for 
each  detection  scheme  and  the  PDF  repre.senting  the  fluctuation  of  the  noise  in  the 
detection  circuit,  tlie  threshold  current  can  be  calculated  for  different  values  of  the 
probability  of  false  alarm. 


5.2.7  The  Probability  of  False  Alarm  for  Incoherent  Detection 

The  equation  for  the  probability  of  false  alarm  for  incoherent  detection  was  given 
in  Eq.  (3.4.4)  as 


Pin.,  =0.5E7?FC| 


^JncoluSN  ^ 


(5.2.1) 


The  dominant  incoherent  noise  when  no  signal  is  present,  w),  is  dark  current 

thermal  noi.se  which  is  1.51x10*^^  A^.  Once  E/Vi/r^  is  chosen  ijcan  be  found.  In  other 
publications  [1,2,4]  probabilities  of  false  alarm  have  been  chosen  between  1()‘*  to  lO'*-^. 
For  this  work,  probabilities  of  false  alarm  of  lO'^,  10’^,  10*^,  10"^,  and  l()‘l^^  were 
cho.sen.  The  complementary  error  function  can  now  be  evaluated  iteratively  until  a 
threshold  current  is  found  that  gives  the  cho.sen  values  for  the  probability  of  fal.se  alarm. 
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The  values  found  for  the  threshold  current  for  the  different  probabilities  of  false  alarm  are 
tabulated  in  Table  5. 


Table  5: 

Threshold  currents  for  different  probabilities  of  false  alarm  for  incoherent  detection. 


Probability  of 
False  Alarm 

if  (nA) 

10-2 

2.860 

10-4 

4.573 

10-^ 

5.846 

10-8 

6.902 

10-10 

7.824 

5.2.2  The  Probability  of  False  Alarm  for  Coherent  Detection 

The  equation  for  the  probability  of  false  alarm  for  coherent  detection  was  given  in 
Eq.  (3.4.6)  as 


(5.2.2) 


Using  the  same  procedure  as  for  incoherent  detection,  the  values  for  the  threshold  current 
were  found  and  are  listed  in  Table  6. 
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Table  6: 

Threshold  currents  for  different  probabilities  of  false  alarm  for  coherent  detection. 


Probability  of 
False  Alarm 

ir  (nA) 

10*-^ 

371.7 

o 

1 

525.7 

10-6 

643.9 

10-^ 

743.5 

10-10 

831.3 

klOLinMMMlMiIOiliTl  jnniTnmnrRTnl 


Now  that  the  threshold  currents  have  been  calculated,  the  probabilities  of  detection 
can  be  calculated  for  each  detection  technique  and  for  each  type  of  target. 


5.3.1  Calculation  of  the  Probability  of  Incoherent  Detection  with  a  Glint  Target 

The  equations  for  calculation  of  the  probability  of  detection  for  incoherent 
detection  with  a  glint  target  were  given  in  Eq.  (3.5.3)  as 


where  /V >0 

(5.3.1) 

where  /r-<'c//m)<0 


where  (/’cin,)  is  the  average  cunent  produced  by  the  signal  incident  on  the  detector,  B  Ls 
the  bandwidth  of  the  detection  circuit,  q  is  the  charge  of  an  electron,  and  {ij^„^sN)  the 

mean  squared  noi.se  current  for  incoherent  detection,  which  is  signal  shot  noise.  These 


56 


wjuations  can  be  used  to  plot  the  probability  of  detection  versus  SNR  fsignal-to-noisc 
ratio)  using  SNR  =  101og|^(/)3,i„,)Y(//„^^j,v)].  the  calculated  mean  squared  dominant  noise 

(Table  4),  the  threshold  currents  for  incoherent  detection  (Table  5),  the  bandwidth  for  the 
incoherent  detection  circuit,  and  by  varying  {ic.,ini)-  This  equation  was  plotted  by  using 

increasing  values  for  which  gave  increasing  values  for  SNR  and  for  the  probability 

of  detection  which  is  easily  seen  in  Fig.  18.  Locking  at  Fig.  18,  as  the  probability  of  false 
alarm  decreases  a  greater  average  SNR  is  needed  to  achieve  the  same  probability  of 
detection  as  for  a  greater  probability  of  false  alarm.  It  can  be  seen  that  for  a  particular 
probability  of  false  alarm,  as  the  average  SNR  increases,  the  probability  of  detection  also 
increases.  It  can  also  be  seen  that  for  a  particular  probability  of  false  alarm  as  the  average 
SNR  increases  the  probability  of  detection  also  increases.  This  makes  sense.  If  the 
average  SNR  increa.ses  then  either  the  noise  has  decreased  (don't  we  wish)  or  the  average 
amount  of  signal  received  has  increased.  If  the  average  amount  of  signal  has  increased 
then  there  would  be  less  of  a  chance  of  missing  a  return,  decreasing  the  probability  of  a 
miss  and  increasing  the  probability  of  detection. 


5.3.2  Calculation  of  the  Probability  of  Incoherent  Detection  with  a  Speckle  Target 

The  equation  to  calculate  the  probability  of  detection  for  tliis  case  was  given  in  Eq. 
(3.5.7)  which  is 


Pd,s  =\\eRFC 


Xr-{{n)+bW) 
^2{{n)-¥bW) 


fcxp(-alV) 

r(M) 


dW. 


(5.3.2) 
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Figure  18:  Curves  representing  the  probability  of  detection  for  incoherent  detection  with 
a  glint  target  for  various  probabilities  of  false  alarm. 


where  Xf  is  the  threshold  number  of  electrons  found  by  dividing  the  threshold  current 
given  in  Table  5  by  ql^t,  where  q  is  the  charge  of  an  electron  and  A?  is  the  transmitted 
pulse  length  (for  this  work  Ar  =  500ns),  M  is  the  number  of  speckle  lobes  received  by  the 
receiver  which  for  equal  entrance  and  exit  aperture  diameters  is  equal  to  4  [15],  (n)  is  the 

mean  number  of  noise  electrons;  that  is,  the  mean  noise  current  for  incoherent  detection 
expressed  in  terms  of  electrons,  W  is  return  energy  incident  on  the  target,  (IV)  is  the 
average  received  energy  from  the  target,  and  fl=M/(w).  Equation  (5.3.2)  was 

integrated  numerically  over  W  and  was  plotted  verses  average  SNR 
(5N/?  =  101og((lV)^/(n)^))  by  using  increasing  values  for  (W)  and  by  using  different 

values  of  the  threshold  current  as  shown  in  Fig.  19.  The  general  trends  of  both  Figs.  IH 
and  19  are  the  same,  but  we  can  compare  the  two  since  incoherent  detection  is  being  u.sed 
in  both  cases  and  the  only  difference  is  the  type  of  target  being  used.  For  a  glint  target  it 
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would  be  expected  that  the  probability  of  detection  will  be  greater  for  a  given  average 
SNR  than  for  a  speckle  target.  This  can  be  seen  in  these  two  plots,  for  a  given  probability 
of  false  alarm,  as  the  average  SNR  increases  the  probability  of  detection  increases  faster 
for  a  glint  target  than  for  the  speckle  target. 
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Figure  19:  Curves  representing  the  probability  of  detection  for  incoherent  detection  with 
a  speckle  target  for  various  probabilities  of  false  alarm. 


5.3.3  Calculation  of  the  Probability  of  Coherent  Detection  with  a  Glint  Target 

An  expression  for  the  probability  of  detection  with  a  glint  target  when  using 
coherent  detection  was  given  in  Eq.  (3.5.9)  as 


(5.3.3) 


Using  the  tlireshold  current,  ij,  given  for  different  probabilities  of  false  alarm  in  Table  6, 
the  dominating  mean  squared  coherent  noise,  (/coa,sw)»  given  in  Table  4  and  by  varying 

iGlint  equation  can  be  plotted  against  the  average  signal-to-noise  ratio  (SNR).  The 
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average  SNR  was  calculated  by  •S'N/?  =  This  equation  for  the 

probability  of  detection  was  plotted  as  described  and  it  is  shown  in  Fig.  20.  The  same 
trends  are  exhibited  in  Fig.  20  as  seen  in  the  previous  plots  for  incoherent  detection. 


Figure  20:  Curves  representing  the  probability  of  detection  for  coherent  detection  with  a 
glint  target  for  various  probabilities  of  false  alarm. 


5.3.4  Calculation  of  the  Probability  of  Coherent  Detection  with  a  Speckle  Target 

The  equation  for  the  probability  of  detection  for  this  case  was  given  in  Eq.  (3.5. 1 1) 
as 


Pdcs  =exp 


{hoH,SN  )  yOiffuu  )  _ 


(5.3.4) 


where  (icoh,sN  )  ^he  mean  squared  noise  current  for  coherent  detection  which  was  given 

in  Table  4  as  1.50x10"^^  A^,  is  the  threshold  current  for  a  particular  probability  of  false 
alarm  and  is  given  in  Table  6,  and  is  the  mean  squared  signal  current  from  the 
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diffuse  target  This  equation  is  easily  plotted  (sec  Fig.  21)  by  varying  )  for  different 
signal  to  noise  ratios  {sNR  =  10  log(^i 

1 

0.98 

I 

I  0,94 
I  0.92 
■5  0,9 
0.88 
0.86 
0.84 
0.82 
0.8 

10  iS  20  25  30  35  40  45  50 

Signal  to  Noise  Ratio  (dB) 

Figure  21:  Curves  representing  the  probability  of  detection  for  coherent  detection  with  a 
speckle  target  for  various  probabilities  of  false  alarm. 


5.3.5  Comparison  of  Coherent  to  Incoherent  Detection  for  a  Glint  Target  Using  the 
Probability  of  Detection 

Now  that  we  have  curves  representing  the  theoretical  probabilities  of  detection  for 
both  coherent  and  incoherent  detection  with  a  glint  target,  they  can  be  compared  to  each 
other  by  plotting  both  Figs.  18  and  20  on  the  same  graph  as  seen  in  Fig.  22.  In  Fig.  22, 
coherent  detection  performs  better  than  incoherent  detection  for  probabilities  of  detection 
greater  that  0.82.  However  the  scale  on  the  horizontal  axis  is  only  one  dB  per  division 
and  the  largest  difference  between  curves  for  a  like  probability  of  false  alarm  is  0.25  dB 
which  is  a  factor  of  1.06.  This  result  was  expected  with  a  glint  target  where  the  return 
signal  is  a  very  strong  deterministic  return.  These  results  suggest  that  for  a  situation  with 
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no  turbulence  the  detection  techniques  should  perform  similarly.  To  understand  the 
surprising  nature  of  this  result  it  is  important  to  remember  that  this  is  for  a  signal  shot 
noise  limited  incoherent  detection  scheme.  For  longer  ranges  this  system  would  no  longer 
be  signal  shot  noise  limited  but  it  would  be  limited  by  the  dark  current  noise  produced  by 
the  dark  current  of  the  photodiode. 
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Figure  22:  Comparison  of  the  probability  of  detection  for  coherent  vs.  incoherent 
detection  with  a  glint  target. 


5.3.6  CompariFon  of  Coherent  to  Incoherent  Detection  for  a  Speckle  Target  Using  the 
Probability  of  Detection 

To  compare  coherent  and  incoherent  detection  for  a  speckle  target  Figs.  19  and  21 
were  combined  and  the  result  is  shown  in  Fig.  23.  As  can  be  seen  in  the  figure,  the 
incoherent  detection  scheme  can  detect  a  target  with  a  lower  signal-to-noise  for  a  given 
probability  of  detection  than  the  coherent  detection  scheme  can.  This  was  an  unexpected 
result.  Since  both  schemes  are  shot  noise  limited  tliey  should  perform  identically.  The 
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difference  in  performance  between  the  two  detection  schemes  is  explained  by  the 
heterodyne  mixing  efficiency  and  speckle  averaging  over  the  receiver  aperture.  In 
coherent  detection  the  mixing  of  the  LO  with  the  received  signal  is  not  perfect  because  of 
distortion  in  the  wave  front  cau.sed  by  propagation  through  the  atmosphere.  This  causes 
only  a  fraction  of  the  return  to  mix  with  the  LO.  Speckle  averaging  occurs  when  the 
receiver  aperture  of  the  system  receives  more  than  one  speckle  lobe.  In  the  incoherent 
case  since  the  detector  dose  not  care  about  the  phase  of  the  return,  all  the  return  energy  is 
seen  by  the  detector.  Although,  this  is  not  the  case  for  coherent  detection.  The  LO  will 
only  mix  with  those  speckle  lobes  that  are  in  phase  with  the  LO.  So  the  incoherent 
scheme  effectively  averages  the  return  speckles  over  the  surface  of  the  detector. 

If  this  system  were  to  be  used  at  larger  ranges  the  system  would  cease  to  be  signal 
shot  noise  limited  and  then  the  system  would  be  thermally  noise  limited.  The  statistics 
would  then  be  different  and  the  coherent  detection  scheme  would  perform  better  than  the 
incoherent  scheme. 
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Figure  23:  Comparison  of  the  probability  of  detection  for  coherent  vs.  incoherent 
detection  with  a  speckle  target. 


CHAPTER  VI 


6.0  Experimfental  Plan 

This  chapter  includes  a  discussion  on  the  data  needed  to  demonstrate  the 
probability  density  functions  used  to  predict  the  probabilities  of  detection,  the  data  needed 
to  verify  the  probability  of  detection,  the  targets  used  to  collect  this  data,  the  shape  of  the 
return  and  how  this  data  was  collected  and  stored. 


6.1  Required  Data 

The  theoretical  predictions  for  the  fluctuation  in  the  detected  output  were  done 
assuming  a  peak  envelope  detected  return  signal  as  will  be  discussed  shortly.  In  order  to 
be  consistent  with  the  theory  the  experimental  return  data  will  also  be  peak  envelope 
detected.  Because  of  the  characteristics  of  the  target  and  the  detection  process  the  peak 
of  the  return  will  fluctuate.  The  distributions  of  these  fluctuations  can  be  seen  by  plotting 
a  histogram  of  the  number  of  pulses  versus  the  amplitude  of  the  peak  of  the  return.  For  a 
large  number  of  shots,  each  bin  represents  the  probability  of  that  amount  of  power 
occurring  as  a  return.  The  distribution  of  the  plotted  data  therefore  represents  the 
probability  function  of  the  data.  By  integrating  the  previously  discussed  PDFs  the 
experimental  data  can  then  be  compared  to  theor>'. 

To  verify  the  probability  of  detection,  data  needs  to  be  collected  for  just  the  noise 
in  the  system  while  the  system  is  running.  Using  this  collected  data  a  threshold  can  be  .set 
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for  a  given  probability  of  false  alarm.  This  threshold  can  then  be  compared  to  the  data 
taken  to  see  what  probability  of  detection  our  system  would  give. 

In  order  to  see  the  fluctuations  just  caused  by  the  target  and  the  detection  scheme 
alone,  the  transmitted  pulse  has  to  be  a  constant  amplitude.  For  this  system,  the  output 
energy  fluctuates  a  good  deal.  An  example  of  the  fluctuation  in  the  output  pulse  can  be 
seen  in  a  histogram  of  the  transmitted  energy  shown  in  Fig.  24.  To  eliminate  this 
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fluctuation  in  the  transmitted  pulse  the  amount  of  energy  transmitted  by  the  laser  will  have 
to  be  recorded  along  with  the  magnitude  of  the  return  pulse.  This  will  allow  for  a 
normalization  of  the  transmitted  pulse  energy  which  will  eliminate  fluctuation  caused  by 
the  fluctuation  in  the  energy  in  the  transmitted  pulse. 
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Figure  24:  Fluctuation  in  output  energy  of  slave  oscillator. 
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The  data  to  be  acquired,  as  just  discussed,  is  the  magnitude  of  the  peak  of  the 
return  and  the  magnitude  of  the  transmitted  pulse  for  both  a  glint  and  a  speckle  target,  as 
detected  by  incoherent  and  coherent  means.  This  section  includes  a  discussion  of  the 
targets  used  to  collect  the  data,  the  format  of  this  data,  and  the  method  used  for  acquiring 
this  data. 

6.2.1  Targets  Used 

Several  different  targets  were  used  throughout  testing.  For  a  speckle/diffuse  target 
a  2'  X  2',  99%  lambertian  target  at  2.09|xm  was  purchased  from  Labsphere  Inc.  Since  this 
target  is  lambertian  it  should  have  nearly  perfect  speckle  target  characteristics.  Along  with 
this  lambertian  target  a  2.5'  x  3.5'  flame  sprayed  Aluminum  target  was  also  used  as  a 
speckle  target.  For  a  glint  target  a  comercube  reflector  was  first  used.  Then  it  was 
realized  that  the  refractive  turbulence  caused  by  temperature  differences  in  the  path  of  the 
beam  were  causing  the  beam  to  wonder  on  and  off  the  comercube.  To  resolve  this 
problem  a  large  retro-reflector  was  needed.  As  a  solution,  65  red  bicycle  reflectors  were 
purchased  and  they  were  mounted  to  a  piece  of  plywood  to  give  a  large  2'  x  2'  glint  target. 

6.2.2  Data  Format 

To  compare  the  data  to  the  theoretical  predictions  the  required  information  is  the 
peak  of  the  return  and  the  peak  of  the  transmitted  energy.  For  coherent  detection  the 
return  waveform,  as  displayed  on  the  Tektronix 's  DSA  602 A  oscilloscope,  output  from  the 
detection  circuit  electronics  is  shown  in  Fig.  25.  The  oscillation  seen  in  the  waveform  is 
near  the  intermediate  frequency  of  27.1  MHz  which  can  be  seen  by  taking  the  Fourier 
Transform  of  the  waveform,  which  is  shown  in  Fig,  26,  The  waveform  used  when  doing 
the  peak  envelope  detection  was  the  waveform  seen  in  Fig.  25.  Normally  an  envelope 
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detector  is  used  to  extract  the  pulse  but,  since  the  system  was  only  operating  at  3  Hz  the 
waveforms  themselves  were  stored  for  later  analysis. 

For  incoherent  detection  the  output  of  the  electronics  is  shown  in  Fig  27.  The 
shape  of  the  return  pulse  is  caused  by  a  differentiation  that  occurs  in  the  detection  circuit 
electronics.  For  this  analysis  the  maximum  of  the  absolute  value  of  the  pulse  was  used  in 
the  data  collection  as  the  peak  of  this  waveform.  Do  to  this,  the  square  of  the  return  was 
taken  and  in  the  data  analysis,  once  the  peak  was  found,  the  .squaie  root  of  that  peak  was 
taken. 


Tlm«  (microseconds) 

Figure  25:  Example  of  a  coherent  return  as  seen  on  the  Tektronix  DSA  6U2A 
oscilloscope. 
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Figure  26:  Example  of  the  Fourier  Transform  of  coherent  detection  return  waveform. 
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Figure  27:  Example  output  of  incoherent  detection  electronics. 

The  final  waveform  that  is  of  interest  is  the  amount  of  transmitted  energy  for  each 
pulse.  This  was  easily  measured  using  the  leakage  from  the  thin  film  Brewster's  angle 
polarizer  shown  in  Fig.  8  using  the  Molectron  JfDlOOO  Joule  meter.  This  meter  has  a  BNC 
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connection  on  the  back  that  is  called  "pulse  output".  This  output  gives  a  voltage  that  is 
proportional  to  the  amount  of  energy  incident  on  the  detector.  The  proportionality 
constant  is  8  V/J  [27],  So  by  dividing  the  peak  of  the  transmitted  pulse  by  8  V/J  the 
average  energy  in  the  transmitted  pulse  will  be  displayed.  An  example  of  the  wavefonn 
stored  from  the  Molectron  is  shown  below  in  Fig.  28.  The  peak  of  Fig.  28  shows  the  peak 
of  the  average  transmitted  energy  for  a  single  pulse  measured  from  the  leakage  through 
the  Brewster's  angle  polarizer.  The  ratio  of  the  transmitted  energy  to  the  leakage  was 
measured  to  be  1.5.3: 1.  So  multiplying  the  peak  of  the  pulse  ui  Fig.  28  by  15.3,  gives 
about  9  mJ  transmitted  energy. 


Figure  28:  Average  Transmitted  pulse  energy  as  displayed  on  the  DSA  602A 
oscilloscope. 
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6,2  J  Data  Collection 

The  best  method  of  data  collection  found  with  the  available  equipment  was  to 
store  each  waveform,  the  transmitted  pulse  energy  and  the  return  signal,  in  tlie  available 
RAM  of  Tektronix's  DSA  602A  oscilloscope.  The  RAM  of  the  oscilloscope  can  hold  up 
to  900  waveforms.  Since  both  the  transmitted  energy  and  the  return  signal  are  being 
stored  the  number  of  actual  data  points  per  data  set  is  only  about  450.  Once  the  memory 
is  full  the  data  can  be  copied  to  a  3.5"  floppy  disk  in  binary  format  for  later  analysis. 


6.3  Experimental  Procedure 

To  take  data  using  the  2.09  |j,m  LADAR  system,  the  system  needs  to  be  up  and 
running.  The  slave  laser  is  the  first  device  to  be  turned  on.  It  is  usually  operated  at  a 
repetition  rate  of  2.8  Hz  with  900  V  applied  to  the  flash  lamps.  It  is  a  good  practice  to 
check  the  alignment  of  this  laser  by  optimizing  the  output  using  the  Molectron  JDIOOO 
energy  meter.  Once  the  output  is  optimized  to  around  15-20  mJ,  Q-switched,  then  it  is 
time  to  turn  on  the  MO,  AOM  and  the  servo  system.  Once  on,  the  servo  system  will  begin 
dithering  the  cavity  trying  to  lock  onto  the  wavelength  of  the  MO.  The  servo  system  may 
also  need  to  be  adjusted  in  order  to  lock  onto  the  correct  mode. 

Once  the  system  is  up  and  running  the  output  beam  can  be  directed  toward  the  test 
field.  The  test  field  used  was  Wright  Field  at  Wright-Patterson  Air  Force  Base.  This 
testing  took  place  fi'om  the  1 1th  floor  of  the  tower  in  building  620,  which  was  built  to  do 
radar  testing  on  Wright  Field.  Using  the  Aerotech  beam  steerer  the  beam  can  be  pointed 
at  the  target  of  interest.  The  targets  of  interest  were  placed  on  the  test  range,  while  other 
targets  are  available  if  desired;  e.g.  tanks,  air  planes,  and  trucks.  Once  the  beam  is 
consistently  striking  the  target  it  is  left  staring  at  the  target  until  both  methods  of  detection 
are  used  and  recorded. 
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All  that  is  left  at  this  point  is  to  get  the  Tektronix,  DSA  602A  oscilloscope  to 
record  the  needed  data.  To  do  this,  the  scope  was  setup  to  trigger  off  of  the  back  scatter 
from  the  secondary  mirror  of  the  telescope  shown  in  Fig.  8,  produced  by  the  transmitted 
pulse.  With  this  oscilloscope  there  is  the  ability  to  display  two  different  wave  forms  on 
two  different  windows.  Therefore  the  return  pulse  is  displayed  in  one  window  and  the 
transmitted  pulse  measured  by  the  Molectron  detector  was  displayed  in  the  other  window. 
Now  that  both  wave  forms  are  displayed,  the  oscilloscope  was  set  to  store  both  wave 
forms  in  RAM.  Once  the  data  was  stored  in  RAM  it  could  be  transferred  to  floppy  disk 
were  analysis  could  be  easily  done  on  the  data  at  a  later  time. 


CHAPTER  Vn 


7.0  ANAI.VSIS  OF  DATA 

This  ciiapter  consists  of  a  discussion  of  peak  detection  of  the  stored  experimental 
data,  manipulation  and  plotting  of  that  experimental  data,  comparison  of  experimental  data 
to  theory  and  determining  the  theoretical  goodness  fit,  explanation  for  the  deviations  from 
theory  and  finally  a  discussion  on  the  probability  of  detection. 

7.1  Peak  Detection 

The  theory  done  to  describe  the  fluctuations  in  the  output  current  from  the 
electronics  assumes  peak  envelope  detection.  As  described  in  the  previous  chapter  the 
entire  wave  form  of  both  the  transmitted  pulse  and  the  return  signal  were  stored.  This 
section  describes  how  the  peak  of  each  return  was  extracted  from  the  stored  data. 

Each  data  set  was  stored  by  the  Tektronix  602A  oscilloscope  on  a  3.5"  disk.  On 
each  disk  are  approximately  900  files  stored  in  sequential  order  using  the  file  names: 
stol.wfb,  sto2.wfb,  sto3.wfb, ...  etc.  Since  the  transmitted  pulse  energy  corresponding  to 
each  return  was  stored,  there  should  be  an  even  number  of  files.  The  file  extension  ".wfb" 
indicates  that  the  file  was  stored  in  binary  format.  Since  both  the  peak  of  the  return  and 
the  peak  of  the  transmitted  energy  are  of  interest  a  program  was  written  in  C  by  Joe 
Gillette,  of  Technology/Scientific  Services,  called  "maxparse,"  that  opened  each  file. 
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searched  for  the  maximum  in  that  file  and  stores  that  value  in  a  user  named  comma 
separated  value  (csv)  file.  Since  the  return  and  the  transmitted  energy  were  stored 
sequentially  alternating  between  the  two,  the  maximums  were  stored  in  two  columns  in  the 
user  named  file. 


7.2  Comparison  of  Data  with  Theoretical  Expectations 

In  this  section  the  methods  used  to  compare  the  experimental  measurements  to  the 
theoretical  expectations  is  discussed.  The  first  section  discusses  hov/  the  transrratted  pui.se 
energy  was  normalized  to  eliminate  the  effects  of  the  fluctuation  in  the  transmitted  energy. 
The  second  section  presents  the  chi^  goodness  of  fit  test  that  will  be  used  to  detemiiije 
whether  the  theory  fits  the  experimental  data.  The  next  four  sections  discuss  the 
probability  distributions  measured  for  each  detection  scheme  for  both  types  of  targets  and 
how  the  theoretical  predictions  fit  those  experimentally  measured  distributions.  Finally, 
the  probability  of  detection  is  discussed. 

7.2.  /  Normalization  of  Transmitted  Pulse 

Once  the  data  set  has  been  reduced  to  a  file  that  is  in  the  co(nma  separated  value 
format  it  can  be  loaded  into  Microsoft's  Excel.  The  rest  of  the  analysis  on  this  data  has 
been  done  in  Excel.  In  order  to  see  the  fluctuation  in  the  detector  current  caused  by  target 
and  detection  characteristics,  the  effects  of  the  fluctuation  in  the  output  pulse  n^«d  to  be 
normalized.  To  do  that  the  measured  voltage  corresponding  to  the  return  is  divided  by  the 
value  recorded  as  the  transmitted  energy.  The  value  recorded  is  not  the  transmitted 
energy  sent  to  the  target  but  it  is  the  leakage  out  of  the  Brewster's  ;uigle  polarizer  which  Ls 
about  1/15.3  of  the  transmitted  pulse.  So  the  actual  transmitted  pulse  is  15.3  times  the 
value  recorded.  This  will  give  some  unu.sual  units  to  the  x-axis  but  that  is  not  of  concern 
since  the  shape  of  the  probability  distribution  is  the  main  concern. 


74 


Now  that  the  pulse  to  pulse  fluctuation  of  the  SO  has  been  accounted  for  the  data 
can  b(5  plotted  in  a  histogram.  The  x-axis  of  the  histogram  is  the  ratio  of  the  received 
signal  voltage  divided  by  the  average  transmitted  energy  and  the  y-axis  is  the  number  of 
times  that  this  values  occurs  during  the  data  set,  which  is  referred  to  as  the  frequency. 

7.2.2  Chi^  Goodness  of  Fit  Test 

When  testing  experimental  data  against  a  hypothesized  probability  density 
function,  the  theoretical  curves  can  be  compared  to  the  experimental  data  and  then  the 
quality  of  the  theoretical  fit  determined.  One  method  of  determining  the  goodness  of  fit  Ls 
a  chi-squared  (x^)  test.  The  statistic  is  defined  as  (see  Appendix  B)  [28] 


where  Oj  is  the  observed  frequency  (experimentally  measured),  Ej  is  the  expected 
(theoretical)  frequency  and  c  is  the  number  of  classes.  Equation  (7.2.1)  repre.sent  a 
random  variable  that  has  an  approximate  x^-distribution  with  c-1  degrees  of  freedom. 
The  number  of  degrees  of  freedom  distinguishes  between  the  members  of  the 
distribution  family.  When  determining  the  number  of  degrees  of  freedom  for  the 
theoretical  distribution  the  equation  is 

Degrees  of  Freedom  =c~\~e  (7.2.2) 

where  e  is  the  number  of  estimated  or  calculated  parameters  in  the  theoretical  distribution 
and  c  in  my  case,  as  will  be  seen  later,  is  the  number  of  bins  in  a  histogram.  Examples  of 
estimated  parameters  would  be  the  mean  and/or  the  variance  used  in  the  theoretical 
distribution  when  they  were  measured  from  the  data. 
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The  distribution  is  a  tabulated  function  that  can  be  found  in  any  statistics  book, 
an  example  of  which  is  shown  in  Appendix  B.  Knowing  the  number  of  degrees  of 
freedom  and  the  level  of  significance,  the  tabulated  value  can  be  compared  to  that 
calculated  for  the  data.  The  level  of  significance  is  the  probability  of  rejecting  the 
hypothesis  even  though  the  hypothesis  is  true.  To  better  explain  the  level  of  .significance 
an  example  is  given  as  follows:  if  the  measured  data  are  samples  of  a  Gaussian  random 
variable,  and  a  decision  threshold  (or  level  of  significance)  of  0. 1  was  selected,  there  is  still 
a  10%  chance  the  data  will  fail  the  test  and  the  hypothesis  will  be  "rejected"  -  that  is,  we 
incorrectly  decide  the  density  is  nol  a  Gaussian.  The  standard  level  of  significance  used 
when  comparing  theoretical  curves  to  experimental  curves  is  10%,  which  has  been  used 
before  in  the  laser  radar  community  and  will  be  used  for  this  work. 

In  order  to  demonstrate  how  the  x^  test  is  carried  out  an  example  follows. 
Assume  that  a  distribution  has  10  degrees  of  freedom  and  a  calculated  x^  value  of  14.6. 
Looking  at  the  tabulated  cumulative  x^-distribution  (Appendix  B)  for  10  degrees  of 
freedom  at  the  0. 1  level  of  significance  (0.90  in  the  table)  the  x^  value  needs  to  be  less 
that  16.0  to  accept  the  theoretical  distribution.  The  calculated  value  was  14.6  which  is 
less  than  16.0  so  this  would  be  accepted  as  a  good  fit. 

7.2.3  Analysis  of  Incoherent  Detection  with  a  Speckle  Target 

In  this  section  the  data  measured  for  incoherent  detection  with  a  speckle  target  will 
be  presented  and  the  process  used  to  fit  the  theoretical  probability  density  function  to  the 
theory  will  be  discussed. 
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Figure  29:  Incoherent  detection  with  the  lambertian  (speckle)  target.  This  data  was 
taken  on  April  12,  1993. 

The  targets  used  in  this  case  were  the  Lambertian  target  and  the  flame  sprayed 
aluminum  target.  Figure  29  is  a  histogram  of  the  returns  from  the  Lambertian  target. 
Before  a  probability  density  function  is  plotted  on  top  of  this  histogram  the  physical 
meaning  of  the  histogram  needs  to  be  understood.  The  height  of  each  bin,  once 
normalized  to  the  total  number  of  returns  measured,  represents  the  probability  of  the 
return  falling  in  that  bin.  The  theoretical  distributions  discussed  earlier  were  probability 
density  functions.  In  order  to  get  a  theoretical  value  for  the  probability,  the  probability 
density  function  has  to  be  integrated  over  the  range  that  represents  the  bin.  Once 
integrated  the  value  remaining  for  each  bin  is  the  probability  of  that  return  occurring. 

Ttie  general  procedure  used  to  fit  the  theory  to  the  data  for  all  four  cases  is  as 
follows: 


1 .  The  mean  and  the  variance  of  the  experimental  datfi  is  calculated. 

2.  A  general  form  of  the  theoretical  distribution  is  found  that  is  a  function  of 
the  mean  and/or  the  variance. 
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3.  The  equations  for  the  mean  and  variance  are  found  in  terms  of  the  variable 
in  the  general  theoretical  distributions. 

4.  The  mean  and  variance  of  the  data  is  used  in  the  equations  for  the  mean 
and  variance  to  solve  simultaneously  for  the  variables  in  the  theoretical 
distributions. 

5.  Using  the  solutions  found  in  step  4  the  general  theoretical  distributions  are 
integrated  over  each  bin  of  the  histogram  and  plotted  versus  the  return. 


The  above  steps  were  used  for  each  detection  scheme  for  each  type  of  target  and  are 
shown  in  the  following  pages. 

The  probability  density  function  which  represents  the  fluctuation  in  the  detector 
current  for  incoherent  detection  with  a  speckle  target  is  given  in  Eq.  (3.3. 10)  as  [15] 


=  J 
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(7.2.2) 


This  equation  is  known  as  a  negative  binomial  distribution.  Fortunately  the  integral  has 
been  simplified  to  a  discrete  form  as  shown  below  [15] 
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(7.2.3) 


where  (yv,)  is  the  average  signal,  M  is  the  number  of  speckle  lobes  seen  by  the  receiver 
and  r  represents  the  gamma  function.  Before  this  equation  can  be  used  to  curve  fit, 
expressions  are  needed  for  the  mean  and  the  variance  of  the  equation.  These  were  not 
given  by  Goodman  [15]  so  they  must  be  derived  as  follows. 

A  general  form  of  a  negative  binomial  distribution  is  [29] 


\w-lj 


(7.2.4) 
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where  the  first  expression  is  the  binomial  coefficient  which  is  defined  as 

'A  xl 

The  mean  is 


w 


0 


(7.2.5) 


(7.2.6) 


and  the  variance  is 


Starting  with  the  binomial  coefficient  in  Eq.  (7.2.3),  Eq.  (7.2.3)  was  put  in  to  the  general 
form  expressed  in  Eq.  (7.2.4).  Through  simple  mathematical  manipulations,  0  was  found 
tobe(yv,)/(Af+(yv,)). 

Knowing  the  mean,  m,  and  the  variance,  Op,  of  the  data,  Eq.  (7.2.6)  and  (7.2.7) 
can  be  solved  for  w  and  0.  The  variance  was  calculated  using 

^1=-^ - ; -  (7.2.8) 

n-l 

where  m  is  the  mean  and  /;}  is  the  n^^  peak  in  a  data  set.  Once  w  and  0  are  known  Eq. 
(7.2,4)  can  be  summed  over  each  bin  of  Fig.  29  to  find  the  theoretical  probability.  A  quick 
way  to  find  the  area  under  a  function  f(x)  between  Xj  and  X2  is  to  perform  calculation  Area 
=  0.5(x2-xi)lf(x2)+f(xj)l.  The  theoretical  probabilities  are  plotted  in  Fig.  30  as  dots  and 
the  curve  plotted  helps  to  visualize  how  a  continuous  probability  function  would  fit  the 
data.  The  chi^  test  was  16.25  for  24  degrees  of  freedom.  This  is  less  than  the  required 
33.2  needed  for  90%  agreement.  This  would  then  be  considered  a  good  fit  to  theory. 
There  are  only  449  data  points  in  this  data  set  which  may  not  be  a  very  good  statistical 
sample.  To  increase  the  number  of  data  points  two  sets  of  data  were  taken  consecutively. 
The  combined  set  of  data  is  shown  in  Fig.  31  along  with  the  theoretical  probability 
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Figure  30:  Theoretical  fit  to  data  taken  using  incoherent  detection  with  the  lambertian 
(speckle)  target  on  April  12,  1993.  There  are  449  data  points  in  this  figure. 

function.  This  fit  has  a  chi^  of  13.6  with  24  degrees  of  freedom.  This  is  in  excellent 
agreement  with  theory.  By  comparison.  Fig.  32  shows  similar  data  taken  from  the  flame 
sprayed  aluminum  target.  The  chi^  value  for  this  case  was  32.6  with  30  degrees  of 
freedom  which  is  a  marginally  good  fit  to  theory.  Another  example  of  data  collected  from 
the  lambertian  target  is  shown  in  Fig.  33.  The  chi^  statistic  was  25.7  with  25  degrees  of 
freedom,  which  gives  a  reasonably  good  fit  to  theory. 


80 


Return  (V)/ Average  Transmitted  Energy  (m  J) 

Figure  31:  Theoretical  fit  to  data  taken  using  incoherent  detection  with  the  lainbertian 
(speckle)  target  on  April  12,  1993.  There  are  898  data  points  in  th-s  figure. 
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Figure  32:  Theoretical  fit  to  data  taken  using  incoherent  detection  with  flame  sprayed 
aluminum  target  on  April  8,  1993.  There  are  898  data  points  in  this  figure. 
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Figure  33:  Theoretical  fit  to  data  taken  using  incoherent  detection  with  the  iambertian 
target  on  April  15,  1993.  There  are  449  data  points  in  this  figure. 


7.2,4  Analysis  of  Incoherent  Detection  with  a  Glint  Target 

For  a  glint  target  a  2"  corner  cube  reflector  was  first  used.  The  data  taken  using 
the  corner  cube  reflector  seem  to  be  weighted  to  the  side  of  a  smaller  return  as  can  be  seen 
in  Fig.  34.  This  behavior  was  noticed  for  a  glint  target  for  both  detection  techniques.  The 
explanation  for  these  results  is  relatively  simple  when  considering  the  circumstances.  The 
temperature  difference  between  room,  temperature  and  the  outside  temperature  was  at 
least  20°F  and  the  porthole  to  the  outside  was  open  allowing  the  air  in  the  building  to  rush 
out  and  mix  with  the  cold  air  right  outside  the  porthole.  This  would  give  a  large  amount 
of  refractive  turbulence  at  the  aperture  of  the  transmitted  beam.  The  effect  of  the 
refractive  turbulence  at  the  transmitting  aperture  would  cause  the  beam  to  wandered  on 
and  off  the  2"  comer  cube  reflector  causing  there  to  be  a  smaller  return  most  of  the  time 
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Figure  34:  Data  taken  using  incoherent  detection  with  the  corner  cube  retlector. 


because  only  part  of  the  beam  was  striking  the  target.  This  would  give  a  large  number  of 
small  returns  as  seen  in  Fig.  34. 

To  correct  this  problem  65  bicycle  reflectors  were  purchased  and  tiled  in  p  2'  by  2' 
square  on  a  piece  of  plywood  and  tins  was  used  to  produce  a  glint  return.  An  example  of 
the  distribution  of  the  fluctuation  in  the  return  with  this  new  target  is  shown  in  Fig.  35. 
This  set  of  data  was  taken  on  April  7, 1993. 

The  theoretical  probability  density  function  for  incoherent  detection  with  a  glint 
target  was  given  in  general  form  in  Eq.  (3.3.2)  as  [2,12,15] 


P,rAx)  = 


ix-mf] 

) 


(7.2.9) 


where  m  is  the  mean  of  the  data  and  is  the  variance  of  the  data.  This  Gaussian  PDF 
can  now  be  integrated  over  each  individual  bin,  which  is  reasonably  easy  because  the 
integral  can  be  reduced  to  an  expression  involving  the  error  function  as  shown  below. 
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Figure  35:  Incoherent  detection  with  the  bicycle  reflector  (glint)  target.  This  data  was 
taken  on  April  7,  1993  at  9:20am.  There  are  449  data  points. 
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This  was  graphed  for  the  mean  and  variance  of  the  experimental  data  and  the  resulting 
curve  can  be  seen  in  Fig.  35. 

This  is  the  only  set  of  data  for  this  case  that  the  theory  actually  fit  One  problem 
that  was  unresolvable  was  that  none  of  the  data  sets  include  more  than  449  data  points. 
There  were  several  attempts  at  combining  consecutive  data  sets  but  the  means  of  each 
dam  set  were  different  enough  that  the  data  sets  could  not  be  combined.  Even  though  a 
larger  "glint"  target  was  being  used  the  same  trend  of  having  an  uneven  distribution 
weighted  to  the  smaller  returns  end  of  the  graph  is  seen  throughout  the  data  sets.  Some 
examples  are  seen  in  Fig.  36  and  37. 

One  observation  about  the  data  is  that  the  distribution  of  the  actual  data  Ls 
narrower  than  the  plotted  theoretical  distribution.  This  is  understandable  since  the  data 


Figure  37:  Incoherent  detection  with  the  bicycle  reflector  (glint)  target.  This  data  was 
taken  on  April  22,  1993. 
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was  taken  in  a  staring  mode,  and  had  everything  been  perfect  there  would  be  a  constant 
return  from  the  glint  target.  But  due  to  turbulence  the  beam  will  wander  some  on  the 
target  giving  some  variance  from  the  mean  return  signal.  If  the  data  was  taken  in  a 
scanning  mode  spatial  averaging  caused  by  the  scanning  would  broaden  the  distribution 
[30].  One  other  possible  explanation  for  the  theory  not  fitting  the  data  is  the  target  used. 
It  is  possible  that  bicycle  reflectors  do  not  give  the  same  statistics  as  does  a  comercube 
reflector  for  glint  target  returns. 


7.2.5  Analysis  of  Coherent  Detection  with  a  Speckle  Target 

Data  was  collected  for  coherent  detection  with  both  the  lambertian  and  the  flame 
sprayed  aluminum  targets.  Both  should  exhibit  diffuse  characteristics.  Shown  in  Figs.  3S 
and  39  are  data  sets  taken  from  the  lambertian  target.  Figure  38  shows  a  histogram  of 
data  taken  on  April  1,  1993.  This  is  a  combination  of  two  data  sets  that  had  the  same 
mean  and  similar  variances  giving  910  data  points. 

The  probability  density  function  representing  tne  fluctuation  in  the  current  in  the 
detection  circuit  electronics  is  a  Rayleigh  PDF  (Eq.  (3.3.12))  as 


This  equation  has  the  general  form  of  [28] 


(7.2.12) 


where  =  (Jp/(2  -f)  and  Op  is  the  variance  of  the  data.  Once  the  variance  of  the  data  is 

known,  Eq.  (7.2.12)  can  be  integrated  from  bin  to  bin  (xi  to  X2)  to  find  the  probability 
function  representing  the  fluctuation  in  the  return.  The  result  of  this  integration  is 
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The  chi^’  st{\tistic  for  this  fit  is  4H.7  with  27  degrees  of  freedom.  This  is  not  an  acceptable 
fit  even  though  it,  at  first  glance,  would  seem  to  be  a  good  fit.  Shown  in  Fig.  39  is  another 
set  of  data  taken  on  April  12,  1993  with  the  lambertian  target.  The  chi^  statistic  for  this 
fit  is  450  which  shows  that  this  is  not  a  good  fit  These  two  figures  are  a  good 
representation  of  all  the  speckle  data  taken. 
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Figure  38:  Coherent  detection  with  the  lambertian  target  taken  on  April  1,  1993.  Shown 
is  the  best  fit  Rayleigh  distribution  with  a  chi^  of  48.7  with  27  degrees  of  freedom. 
(Needs  to  be  less  than  36.7). 

One  explanation  for  the  deviation  from  theory  that  has  already  been  mentioned  in 
the  previous  section  is  the  effects  of  turbulence.  With  the  temperature  difference  between 
outside  and  inside  the  tower  the  turbulence  right  outside  the  porthole  has  a  large  effect. 
Some  work  has  been  done  by  Jeffrey  Shapiro,  Sun  T.  Lau  and  David  M.  Papurt  [30,31]  to 
describe  the  effects  of  turbulence  on  target  returns  and  on  target  statistics.  In  their  work 


87 


f31]  they  give  the  probability  density  functions  for  both  a  glint  and  a  speckle  target 
including  the  effects  of  turbulence  using  log-normal  distributions  which  are  functions  of  a 
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Figure  39:  Coherent  detection  with  the  lambertian  target  taken  on  April  12,  1993.  There 
are  898  data  points.  This  is  the  best  fit  Rryleigh  probability  distribution  which  gives  a 
chi2=450. 

which  is  used  as  a  turbulence  factor  that  ranges  from  0  (no  turbulence)  to  0.5  (saturated 
scintillation).  In  this  work  they  also  give  the  cumulative  probability  distribution  function 
for  each  case.  The  cumulative  probability  distribution  for  a  speckle  target  is  given  as 

P,  (X)  =  1  -  Fr[x  exp(4a'  ),0; a];  X  >  0  (7.2.1 4) 

where  the  log-normal-density  frustration  function,  FrlJ,  is  given  as 

(  ^ 

Fr(a,0;c)=  f  dr(27tc^)~^exp - 5 —  exp(-ae‘')  (7.2.15) 

2c‘ 

and 
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where  i  is  the  return  signal  current  and  is  the  mean  of  the  square  of  the  return  signal. 

Since  the  cumulative  probability  distribution  is  the  integration  from  negative 
infinity  to  x  the  equation  can  be  used  to  find  the  probability  of  that  return  occurring  for 
each  bin.  The  calculation  of  the  frustration  function  was  performed  using  a  macro  in 
Microsoft  Excel.  This  analysis  was  done  for  the  collected  data  and  the  difference  this 
makes  can  be  seen  in  Figs.  40  and  41  which  are  the  same  data  sets  as  seen  in  Figs.  38  and 
39  but  with  different  theoretical  fits.  In  Fig.  40  the  chi^  statistic  was  22.5  with  27  degrees 
of  freedom  which  gives  very  good  agreement  with  theory. 


Figure  40:  Data  taken  April  1,  1993  with  log-normal  distribution  used  as  best  fit  curve,  a 
=  0.05,  chi2=22.5  with  27  degrees  of  freedom, 


Figure  41  shows  the  data  taken  on  April  12  with  the  best  fit  log-normal  distribution.  The 
chi^  statistic  was  19.2  with  29  degrees  of  freedom  which  is  excellent  agreement  with 
theory. 


7.2.6  Analysis  of  Coherent  Detection  with  a  Glint  Target 

For  this  case  the  data  was  collected  from  the  target  composed  of  the  bicycle 
reflectors.  The  probability  density  function  representing  the  fluctuation  in  the  output 
current  of  the  detector  electronics  is  represented  by  a  Rician  PDF.  A  general  form  of  the 
Rician  density  is  given  as  [2] 
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Figure  41:  Data  taken  April  12,  1993,  with  best  fit  log-normal,  a=0.055,  chi2=19.2  with 
29  degrees  of  freedom. 


where  is  the  variance,  which  is  related  to  the  variance  of  the  data  represented 

by  [32] 


=  2ct  -f- 
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and  the  mean  is 
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wherej?is/o,i„,/0f. 


To  plot  the  theoretical  curves  with  the  collected  data  the  mean  and  the  variance  of 
the  data  will  be  used  in  conjunction  with  Eqs.  (7.2.16)  and  (7.2.17)  to  give  two  equations 
with  two  unknowns  (i.e.,  iQlint  ^)*  Once  the  values  for  these  variables  are  found  they 
can  be  substituted  into  Eq.  (7.2.15)  and  this  equation  can  be  integrated  over  each  bin  to 
give  the  probability  of  that  return  occurring.  For  the  data  collected,  complex  numbers 
were  calculated  in  the  process  of  trying  to  solve  Eqs.  (7.2.16)  and  Eq.  (7.2.17),  which 
suggests  that  the  Rician  distribution  will  not  fit  the  data.  Some  examples  of  the  data 
collected  are  shown  in  Figs.  42-44.  The  histograms  of  the  data  look  similar  which 
suggests  that  the  same  problem  exists  for  aU  the  data  sets.  An  explanation  for  the 
deviation  from  the  theory  is  refractive  turbulence. 

In  the  previous  section  refractive  turbulence  was  shown  to  have  an  effect  on  the 
probability  distribution  of  the  return  data.  In  Shapiro  and  Lau's  paper  [29]  they  give  an 
equation  that  represents  the  cumulative  probability  function  for  the  fluctuation  of  the 
return  signal  for  coherent  glint  target  detection  with  atmospheric  turbulence  which  is 
expressed  as 


(7.2.18) 


where 


(7.2.19) 


(7.2.20) 


NormalizecI  Frequefx:y  (Probability)  Normalized  Frequency  (Probability) 
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0.08 


0.00  0.06  0.23  0.61  0.90  1.41  2.03  2.77  3.61  4.57  5.64  6.83 

Return  (V)/ Average  Transmitted  Energy  (m  J) 

Figure  42:  Data  taken  with  a  glint  target  on  April  S,  1993. 


0.12 

0.1 


0.00  0.10  0.42  0.94  1.67  2.61  3.76  5.12  6.69  8.46  10.45 

Return  (V)/ Average  Transmitted  Energy  (m  J) 

Figure  43:  Data  taken  with  a  glint  target  on  April  7,  1993. 


0.00  0.15  0  60  1.36  2.40  3.74  5.39  7.34  9.59  12.13  14.98 

Return  (V)/ Average  Transmitted  Energy  (mJ) 

Figure  44:  Data  taken  with  a  glint  target  on  April  8,  1993.  There  are  1353  data  points  in 
this  histogram. 

i  is  the  return  signal  current,  is  the  mean  square  return  signal  and  a  is  the  turbulence 

factor.  Since  Eq.  (7.2.18)  is  a  cumulative  distribution  it  can  be  used  to  find  the  theoretical 
value  of  the  probability  over  each  bin.  This  looks  and  sounds  complicated  but  it  is 
relatively  easy  since  the  Q(t)  function  is  in  the  form  of  the  error  function  (ERF).  The 
expression  used  to  find  the  probability  over  each  bin  (with  boundaries  X]  and  X2)  is  shown 
below. 


(7.2.21) 
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Figure  45:  Theoretical  probability  distribution  for  coherent  detection  with  a  glint  target 
for  data  taken  on  April  5.  Chi2=37.0  with  30  degrees  of  freedom;  good  fit  to  data. 


Return  (V)/ Average  Transmitted  Energy  (m  J) 

Figure  46:  Theoretical  probability  distribution  for  coherent  detection  with  a  glint  target 
for  data  taken  on  April  7.  Chi2=41.8  with  29  degrees  of  freedom;  this  is  a  poor  fit  to 
theory. 
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Equation  (7.2.21)  was  plotted  for  the  data  sets  shown  in  Figs.  42-44.  Shown  in 
Fig.  45  is  the  data  taken  on  April  5^^  with  equation  (7.2.21)  plotted  with  a  a=0.0992 
which  gives  a  chi^  of  37.0  with  30  degrees  of  freedom.  This  is  a  good  fit  to  the  theory.  In 
Fig.  46,  the  best  fit  curve  gave  a  o=0. 1,  chi^  =  41.8  with  29  degrees  of  freedom.  This  is  a 
poor  fit  to  the  theory.  Figure  47  is  a  combination  of  three  data  sets  (1353  data  points) 
taken  on  April  8th.  The  turbulence  factor  for  the  best  fit  was  found  to  be  0.104,  which 
gave  a  chi^  =  76  with  49  degrees  of  freedom.  This  is  also  a  poor  fit  to  theory. 

These  graphs  demonstrate  that  turbulence  is  a  factor  that  is  contributing  to  the 
shape  of  the  probability  curves.  The  theory  that  Shapiro  and  Lau  demonstrated  was 
demonstrated  for  a  10.6  |im  lidar  system.  Since  we  are  working  at  2.09  pm  the  turbulence 
effects  are  5  times  greater.  Another  factor  that  may  effect  the  shape  of  the  distribution  is 
if  the  bicycle  reflectors  do  not  make  a  good  glint  target. 


Figure  47:  Coherent  detection  data  taken  from  a  glint  target  on  April  8^^.  The 
theoretical  fit  gives  a  chi2=76  with  49  degrees  of  freedom.  This  is  a  poor  fit  to  the 
experimental  data. 
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7.2.y  Analysis  of  the  Probability  of  Detection 

To  verify  the  probability  of  detection,  data  was  taken  of  just  the  noise  present  in 
the  system.  The  data  sets  contained  898  data  points  each  and  there  were  6  different  sets 
of  data  taken.  With  each  noise  data  set  a  threshold  was  chosen  to  give  a  particular 
probability  of  false  alarm.  For  example,  one  sample  of  noise  may  have  values  of  2,  5, 7,  4, 
and  9.  If  a  threshold  of  8  is  chosen  the  probability  of  false  alarm  would  be  0.2  (1/5)  since 
there  is  one  value  of  noise  greater  than  the  dtfeshold  for  the  5  pieces  of  data. 

Now  that  a  threshold  is  set  the  data  taken  is  compared  to  that  threshold.  The 
number  of  data  points  greater  than  the  threshold  are  counted  and  that  number  is  divided  by 
the  number  of  data  points  taken.  This  will  give  the  probability  of  detection.  Since  the 
probability  of  false  alarm  is  known  the  amount  of  noise  and  signal  is  also  know  so  the 
average  signal-to-noise  ratio  can  be  calculated  and  then  the  calculated  probability  of 
detection  can  be  compared  to  the  theory  for  that  false  alarm  rate  and  that  SNR.  Because 
of  time  constraints  on  use  of  the  experimental  equipment  and  ranges,  not  enough  data  was 
taken  to  verify  the  probability  of  detection  for  the  glint  targets,  but  sufficient  data  was 
taken  to  verify  the  probabilities  of  detection  for  a  speckle  target. 

For  coherent  detection  with  a  speckle  target  the  results  of  the  measurements  taken 
are  shown  in  Table  7  for  a  probability  of  false  alarm  of  10'2, 

Table  7: 


Comparison  of  calculated  and  measured  probabilities  of  detection  for  coherent  detection 

with  a  speckle  target. 


RIffli 

Probability  of  Detection  Measured 

Probability  of  Detection  Predicted 

30.33 

0.9978 

0.9958 

30.63 

0.9955 

0.996 

29.68 

0.9978 

0.995 
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For  incoherent  detection  with  a  speckle  target  the  results  of  the  measurements 
taken  are  shown  in  Table  8  also  for  a  probability  of  false  alarm  of  10‘2. 

Table  8: 


Comparison  of  calculated  and  measured  probabilities  of  detccdon  for  incoherent  detection 

with  a  speckle  target. 


HSHi 

Probability  of  Detection  Measured 

Probability  of  Detection  Predicted 

18.81 

0.9978 

0.970 

21.41 

0.9911 

0.989 

20.53 

0.9978 

0.985 

22.23 

0.9978 

0.992 

21.75 

0.9978 

0.991 

The  limited  amount  of  data  available  for  both  coherent  and  incoherent  detection 
schemes  seem  to  demonstrate  a  rather  close  comparison  between  the  theory  and  the  data. 


/ 

1. 


CHAPTER  Vni 


g.P  CONCLUSIONS  ANP  RECOMMENDATIONS 

The  basis  used  to  compare  coherent  to  incoherent  detection  for  a  solid  state  2,09 
|xm  LADAR  system  was  the  probability  of  detection.  In  order  to  find  the  probability  of 
detection  the  dominating  noises  in  the  LADAR  system  had  to  be  determined.  It  was 
determined  that  signal  shot  noise  was  the  primary  noise  source  for  incoherent  detection 
and  for  coherent  detection  the  dominating  noise  was  found  to  be  thermal  noise.  Since  the 
desire  for  an  optimized  coherent  detection  scheme  is  to  have  the  system  shot  noise  limited, 
the  amount  of  LO  power  incident  on  the  detector  was  increased  until  the  shot  noise 
dominated  the  thermal  noise  in  the  system.  Knowing  the  dominating  noise  source  for  each 
detection  scheme  the  threshold  currents  for  selected  probabilities  of  false  alarm  were 
determined.  Using  the  threshold  currents  calculated,  the  probabilities  of  detection  for  each 
detection  scheme  were  plotted  for  different  probabilities  of  false  alarm. 

To  calculate  the  probabilities  of  detection  for  both  detection  schemes  for  both  glint 
and  speckle  targets  the  corresponding  probability  density  functions  (PDF's)  representing 
the  fluctuation  in  the  current  from  the  detector  were  identified.  To  verify  that  these  PDFs 
were  accurate,  data  was  collected  of  the  fluctuation  in  the  detection  circuit  electronics 
fi'om  pulse  to  pulse.  An  attempt  was  then  made  to  fit  the  predicted  theory  to  the 
experimental  results. 

For  incoherent  detection  with  a  speckle  target  the  theoretical  probability 
distribution  used  to  find  the  probability  of  detection  fit  well  with  the  experimental  data. 
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This  confirmed  that  the  negative  binomial  distribution  was  the  correct  distribution  to 
represent  the  statistics  of  the  return  for  a  speckle  target  for  incoherent  detection. 

When  using  a  glint  target  with  incoherent  detection  the  results  were  not  as  good. 
There  was  one  data  set  that  fit  the  theory  according  to  the  chi^  statistic.  The  rest  of  the 
data  gave  distributions  that  were  narrower  than  the  theoretical  distributions  that  were 
supposed  to  fit  that  data.  An  explanation  for  this  discrepancy  involves  the  method  used  to 
collect  the  data.  The  data  was  collected  in  a  starring  mode  so  if  there  were  no  other 
effects  to  degrade  the  return  (turbulence,  vibrations,  etc.)  the  return  would  be  a  constant 
value  from  pulse  to  pulse.  Since  the  world  is  not  perfect  those  effects  are  present  so  there 
is  variation  in  the  return,  but  the  return  is  not  entirely  random  because  the  beam  is  still 
striking  the  same  general  location  on  the  target.  This  is  seen  by  looking  at  the  data 
collected  and  noticing  that  the  distribution  of  the  experimental  data  is  narrower  than  the 
theoretical  distribution.  One  recommendation  to  solve  this  problem  is  to  scan  the  beam 
over  an  area  that  includes  the  target  and  then  take  the  value  that  represents  the  largest 
return  from  that  fi'ame.  This  would  give  a  more  random  sample. 

For  coherent  detection  with  a  speckle  target  the  data  did  not  fit  the  Rayleigh 
distribution  as  predicted  by  the  theory.  Under  the  conditions  which  data  was  taken  the 
suspected  problem  is  turbulence.  The  temperature  difference  between  outside  and  inside 
the  test  tower  was  at  least  20°F.  The  porthole  used  to  direct  the  beam  outside  was  not 
covered  so  the  air  in  the  building  was  rushing  out  of  the  room  and  mixing  with  the  outside 
air  immediately  outside.  This  would  give  a  good  deal  of  refractive  turbulence  right 
outside  the  porthole.  Turbulence  at  the  transmit  aperture  is  the  worst  place  to  have 
turbulence  in  a  radar  situation.  Work  done  by  J.H.  Shapiro,  and  S.T.  Lau.  [31]  involved 
describing  the  effects  of  turbulence  on  the  statistics  of  target  returns,  'Fhey  gave  two 
probability  density  functions  that  represent  the  fluctuation  in  the  return  for  a  glint  and  a 
speckle  target  which  included  a  term  representing  the  turbulence  in  the  system.  The 
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equation  for  a  speckle  target  was  plotted  using  the  data  collected  and  there  was  a  verj/ 
good  fit  to  the  theory. 

For  a  glint  target  when  using  coherent  detection  the  fluctuation  in  the  return 
should  be  a  Rician  distribution.  In  order  to  plot  the  theoretical  curve  to  the  experimental 
data  two  equations  with  two  unknowns  have  to  be  solved  to  find  the  parameters  to  graph 
the  equation.  With  the  collected  data  the  result  was  always  complex  which  showed  that 
the  theory  did  not  fit  to  the  data.  Again  it  was  assumed  that  the  reason  was  turbulence. 
Using  the  probability  density  function  given  by  J.H.  Shapiro,  and  S.T.  Lau.  [31]  which 
included  terms  representing  the  turbulence,  a  couple  sets  of  data  were  found  to  fit  to  that 
theory,  but  not  very  well.  Shapiro  and  Lau  demonstrated  for  a  staring  system  that  the 
statistics  would  not  fit  the  theoretical  distribution  and  to  get  the  theoretical  statistics  a 
scanning  system  is  required. 

Using  the  LADAR  system  a  couple  points  on  the  probability  of  detection  curves 
for  speckle  targets  for  coherent  and  incoherent  detection  were  verified.  Because  of  the 
lack  of  time  data  for  the  glint  target  case  is  not  available. 

A  few  suggestions  for  future  work  are:  first  of  all,  data  needs  to  be  taken  in  a 
scanning  mode  which  would  allow  a  more  random  sample  of  target  returns  and  this  would 
give  a  better  representation  of  the  statistics  of  interest.  A  second  suggestion  is  to  take  the 
data  with  something,  a  piece  of  glass  or  plastic,  over  the  porthole  decreasing  the  effects  of 
turbulence  at  the  transmitting  aperture  or  even  opening  the  entire  room  to  the  outside. 
The  third  suggestion  is  to  work  with  the  statistics  involved  for  incoherent  detection  when 
the  system  is  not  signal  shot  noise  limited  and  see  how  the  detection  techniques  compare. 
The  fourth  suggestion  is  to  have  a  better  method  of  collecting  the  data.  The  method  used 
was  the  best  available  with  the  time  and  equipment  available.  Before  the  current  system 
can  be  used  for  further  work,  some  experimental  work  needs  to  be  done  to  determine  why 
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»*’*’■  system  is  losing  efficiency  and  not  transmitting  as  much  energy  as  before.  ITiere  also 
s  to  be  improvements  in  the  consistency  of  the  injection  seeding  of  the  slave  cavity. 

8.1  LIST  OF  ORIGINAL  CONTRIBUTIONS 

•  First  to  investigate  detection  statistics  with  a  coherent  solid  state  laser  radar 
system. 

•  First  to  investigate  detection  statistics  at  2  pm. 

•  First  to  compare  coherent  detection  to  incoherent  detection  at  2  pm 


APPENDIX  A 


An  Understanding  of  Shot  Noise  Using  Poisson  Statistics 

Shot  noise  is  caused  by  fluctuation  in  the  current  due  to  the  discrete  nature  in 
which  charge  carriers  are  produced.  To  discuss  this  fluctuation  there  are  a  few 
assumptions  that  have  to  be  made.  With  those  assumptions  a  brief  description  of  the 
process  that  occurs  during  detection  is  given  and  it  is  described  using  a  Poisson  probability 
distribution.  It  is  then  shown  that  the  mean  squared  fluctuation  of  a  Poisson  probability 
distribution  is  equal  to  the  mean.  This  result  is  related  to  current  from  our  ideal  detector 
and  the  result  is  an  expression  for  the  mean  squared  fluctuation  in  the  current  due  to  the 
detection  process  alone. 

An  ideal  photon  detector  is  assumed.  This  ideal  photon  detector  has  no  leakage 
current  when  there  is  no  light  incident  upon  it.  The  beam  incident  on  the  detector  is 
assumed  noiseless,  single  frequency,  and  of  a  constant  power.  This  would  lead  to  the 
thought  that  the  output  current  would  also  be  constant,  though,  this  is  not  the  case.  The 
depth  at  which  the  photons  travel  into  the  detector  before  they  are  absorbed  and  produce 
the  emission  of  an  electron  is  discrete  and  random  which  means  that  the  emission  of  the 
electrons  is  also  a  random  discrete  process.  A  random  discrete  process  can  be  described 
by  a  Poisson  probability  distribution  as  [32] 

=  (A.I) 

n\ 
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where  P(n)  is  the  probability  that  n  photons  will  be  detected  in  a  sample  time  t,  and  <n> 
is  the  average  number  of  photon  events  oc(;urring  in  time  x. 

The  definition  of  noise  is  the  mean  squared  fluctuation,  or  the  variance.  The  variance  can 
be  calculated  by  [12] 

o' =((«-(«»"),  (A.2) 

which  can  be  simplified  to[32] 

o’  =((n-(B))’^  =  (n’-2n(n)+(n>’) 

=  {«’)-<2B<n»+{(«)’) 

=  (n’)-(2(n>’)+((n)’) 

=  («’>-<«>’ 


This  result  can  now  be  used  to  show  that  the  variance  of  a  Poisson  distribution  is 
equal  to  the  mean.  To  begin,  a  simple  mathematical  identity  is  used. 

=n+n{n-i)  (A.4) 

The  mean  of  Eq  (A.4)  is 

(n  )  =  (n+«(n-l)) 

'  '  ,  \  (A.5) 

=  {n)+{n{n-l)). 

For  the  second  term  of  that  result,  the  following  statistical  definition  can  be  substituted 
[32] 

{n{n  - 1))  =  X  n(/i  -  \)P{n). 

wO 


(A.6) 
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Substituting  Eq.  (A.6)  into  Eq.  (A.5)  we  obtain 

=  (A.7) 

Eq.  (A.7)  can  now  be  simplified  by  noticing  that  n(n-J)  is  the  product  of  the  last  two  terms 
in  the  n!  series,  and  that  the  summation  can  begin  with  n=2  since  the  n=l  and  n=0  terms 
are  both  zero.  In  addition  the  constant  factor  <n>^  can  also  be  pulled  out  of  the 
summation  to  give 

=  (A.8) 

We  now  make  a  change  of  variable  to  m=n-2  and  reduce  Eq.  (A.8)  as  follows  [32] 

„m0 

=  ('■)+('■>' 2 (A.9) 

m*0 

=  (n>+(n)^ 


Where  we  recall  that  the  summation,  over  all  m,  of  the  Poisson  distribution  is  exactly  one. 
This  result  can  then  be  substituted  into  Eq.  (A.3),  which  gives 


a*=(«> 


(A.10) 


Thus  the  variance  of  a  Poisson  distribution  is  equal  to  the  mean.  This  result  will  now  be 
used  to  And  the  mean  squared  fluctuation  in  the  detector  current. 
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If  the  current  is  monitored  at  the  output  of  the  detector  for  a  period  of  time,  t  the 
current  measured  would  be 


(A.ll) 


where  n  is  the  number  of  electrons  generated  by  the  detector  in  time  t  and  q  is  the  charge 
on  an  electron.  The  average  current  is 


(A.  12) 


The  mean  squared  fluctuation  in  the  current,  which  can  also  be  called  the  mean  squared 
noise  current,  is 


(A.  13) 


The  second  term  in  Eq.  (A.  13)  is  just  the  variance,  which  has  been  shown  above  is  the 
mean.  Using  Eq.  (A.  12)  the  mean  can  be  substituted  and  the  mean  squared  noise  current 
becomes 

=  (A.14) 

where  I  is  used  as  the  average  current  Hnally  the  sample  time  x  is  equal  to  one  over  two 
times  the  bandwidth.  This  can  be  substituted  into  Eq.  (A.14)  to  give  the  final  form  of  the 
mean  squared  noise  current  commonly  known  as  shot  noise  [12];  that  is 
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(4)  =  2<7/B 


(A.  15) 
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APPENDIX  B 

Chi2  Goodness  of  Fit 

A  goodness  of  fit  test  determines  whether  a  set  of  data  can  be  viewed  as  a  random 
variable  with  a  given  distribution.  To  test  whether  measured  distribution  actually  fits  the 
proposed  theoretical  distribution  a  level  of  significance  (a)  is  chosen.  For  this  work  the 
level  of  significance  used  is  0.1,  which  means  there  is  a  (l^allOO  =  90  percent  confidence 
that  the  measured  distribution  fits  the  proposed  theoretical  distribution.  One  method  of 
determining  the  goodness  of  fit  is  a  chi-squared  (x^)  test.  The  statistic  is  defined  as 
[28] 


where  Oj  is  the  observed  firequency  (experimentally  measured),  Ej  is  the  expected 
(theoredcal)  frequency  and  c  is  the  number  of  classes.  Equation  (B.l)  represents  a 
random  variable  that  has  an  approximate  x^-distribution  with  c-1  degrees  of  freedom. 
The  number  of  degrees  of  freedom  distinguishes  between  the  members  of  the  x^- 
distribution  family.  When  determining  the  number  of  degrees  of  freedom  for  the 
theoretical  distribution  the  equation  is 

Degrees  of  Freedom  =c-l-e  (B.2) 

where  e  is  the  number  of  estimated  or  calculated  parameters  in  tlie  theoretical  distribution 
and  c  in  my  case,  as  will  be  seen  later,  is  the  number  of  bins  in  a  histogram.  Examples  of 
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estimated  parameters  would  be  the  mean  and/or  the  variance  used  in  the  theoretical 


distribution  when  they  were  measured  firom  the  data. 

The  distribution  is  a  tabulated  function  that  can  seen  listed  below  dependent  on 
the  number  of  degrees  of  freedom  and  the  level  of  significance.  To  demonstrate  how  to 
use  this  table  we  will  use  a  0.1  level  of  significance  and  24  degrees  of  freedom.  Looking  at 
Table  B1  the  number  of  degrees  of  freedom  (D  of  F)  is  shown  in  the  left  most  column  and 


1-Oevel  of  significance)  is  shown  in  the  top  row.  To  use  this  table  follow  down  the  D  of  F 


column  until  24  degrees  of  freedom  is  found,  then  go  across  to  the  column  that  is  0,9, 
which  gives  a  of  33.2.  So  in  order  to  accept  a  theoretical  fit  to  a  set  of  data  the  of 
that  data  must  not  be  more  that  33.2. 
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